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ABSTRACT 


Software for data archiving and data display was developed for use on a 
Digital Equipment Corporation (DEC) PDP-11 /34A minicomputer for use with the 
JPL-designed flux mapper. The flux mapper is a two-dimensional, high radiant 
energy scanning device designed to measure radiant flux energies expected at 
the focal point of solar parabolic dish concentrators. Interfacing to the DEC 
equipment was accomplished by standard RS-232C serial lines. The design of 
the software was dictated by design constraints of the flux-mapper controller. 
Early attempts at data acquisition from the flux-mapper controller were not 
without difficulty. Time and personnel limitations resulted in an alternative 
method of data recording at the test site with subsequent analysis accomplished 
at a data evaluation location at some later time. Software for plotting was 
also written to better visualize the flux patterns. Recommendations for 
future or alternative development are discussed. A listing of the programs 
used in the analysis is included in an appendix. 
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SECTION I 


INTRODUCTION 


A. OVERALL DESIGN PHILOSOPHY 

The flux mapper is a three-dimensional scanning system to measure the 
high radiant flux levels expected at the focal point of a solar parabolic dish 
system. The scanning, measurement, and initial storage of the data are 
handled by the flux-mapper controller. Software was written to enable a 
Digital Equipment Corporation (DEC) PDP-11/34A minicomputer to archive the 
data for long-term storage and to display the data collected from the 
flux-mapper controller in different formats. 

The overriding system requirement for the software was compatibility 
with the output of the flux-mapper controller. The design, fabrication, and 
software control of the flux mapper were carried out by JPL personnel. The 
flux-mapper controller is a microprocessor-based system using read-only memory 
(ROM) to store the acquisition and output routines. Consequently, changes in 
the software routines from the flux-mapper controller are more difficult to 
achieve. It was decided that the requirement for compatibility would rest 
with the PDP-11/34A because the programming was to be done using standard 
Fortran IV language and compiler. Therefore, changes in gathering of the data 
from the controller could be made quickly. 

Another important design requirement was ease of use in the field. The 
flux mapper was to be used at the Parabolic Dish Test Site (PDTS) located at 
the JPL Edwards Test Station (ETS) in the high desert, 120 km northeast of JPL 
(in Pasadena). The personnel at the PDTS were not extensively trained in 
minicomputer programming. Also, the small size of the staff precluded the 
availability of an individual dedicated to the minicomputer system. 

Therefore, the staff at the PDTS did not have the time nor the expertise 
to perform significant minicomputer tasks. 

Finally, it was felt that error checking of the incoming data from the 
flux-mapper controller should be included. Incorrect characters would 
obviously be detrimental to output displays and listings as well as to the 
archived data. Because the PDP-11 /34A minicomputer would receive data from 
the flux-mapper controller, error checking could only be done as the data were 
received. This precludes the use of standard methods such as checksums; 
therefore, a different approach had to be used. 

The three above-mentioned overall design requirements were the basis by 
which all software was designed. These initial decisions were made during the 
end of 1979 when discussions were first taking place regarding the flux mapper 
and its interface with the data-acquis it ion minicomputer system at the PDTS. 
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B. 


DATA- ACQUISITION HARDWARE 


The data-acquisition hardware available at the PDTS consisted of a 
Digital Equipment Corporation PDP-11/34A minicomputer with two removable RK05 
disk drives and 256 kilobytes of internal memory. The system also had as 
supporting peripherals a Kennedy Model 9100 magnetic tape drive, a Versatec 
Model 1100 printer /plotter , and a Control Data Corporation Model 9766 
removable disk drive. In addition, there were interfaces for the various 
terminals used as well as additional serial ports for access by other devices 
such as the data loggers . The entire minicomputer system was housed in a 
mobile trailer located adjacent to the PDTS control room at the test site. 

All connections between the minicomputer, the terminals, and data loggers 
(with the exception of the console terminal) were by RS-232C standard serial 
interfaces. The console terminal communicated by the standard 20-milliamp 
current loop. 


C. GENERAL PRACTICAL RESTRICTIONS 

A significant restriction that was not expected occurred as a result of 
the operating system used. The multi-user system supplied by Digital 
Equipment Corporation (DEC) called RSX-11M was used during the entire period 
of development. Unlike the DEC single-user system, the RT-11, their 
multi-user system would not permit the use of a ring buffer to store input 
data. During each cycle, the operating system polls all input devices to 
determine if a task has been initiated. During this polling procedure, inputs 
to other ports are not placed into ring buffers for temporary storage until 
they are polled in turn. This resulted in the possibility of losing input 
characters from a device while the operating system was polling other devices. 

With this problem in mind, it became clear that during transfer of data 
from the flux-mapper controller to the PDP-11/34A the possibility of data loss 
would be great if this polling mechanism were left intact. This could cause a 
significant restriction of implementation of the data-acquisition and 
archiving software. 
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SECTION II 


INITIAL DATA-LOGGING ATTEMPTS 


A. DESIGN PHILOSOPHY 

During the initial phases of the development of the necessary software, 
the same general design criteria stated previously were used. Compatibility, 
ease of use, and error checking were of foremost importance in the design of 
the software. 

During software development, accessibility of the flux mapper for 
testing was limited. Therefore, much of the development of the software had 
to be done without access to the actual hardware. This necessitated the 
inclusion of ease of hardware interfacing as a design criterion. And finally, 
because the equipment was to be used at the PDTS, where the personnel did not 
have minicomputer or microprocessor expertise, minimal interfacing in terms of 
training and operator instructions was a desirable criterion. 


B. AVAILABLE HARDWARE 

The available hardware was discussed previously. The flux-mapper 
controller output, basically an RS-232C serial interface line, was connected 
to one of the serial inputs of the minicomputer. 


C . SOFTWARE WRITTEN 

The initial software written consisted of a data-acquisit ion program and 
a data-output program. Both were written in Fortran IV for the DEC PDP-11/34A 
minicomputer using the RSX-11M operating system. 

The data-acquisition software was written to read the seven-bit ASCII 
characters that are transmitted from the flux-mapper controller. The format 
for the data from the flux-mapper controller consists of three types of 
records: a header record, a set of data records, and an end record. 

Each record was identified by a one-character ASCII identifier, followed 
by the data in a fixed structure. The structures of each type of record are 
given in Table 1. Note that there was only one header record and only one end 
record, but there could be many data records. 

Each data record represented one traverse of the radiometer probe across 
the flux mapper. At the end of the traverse, the probe incremented in the 
perpendicular direction by a preset amount and continued its traverse in the 
opposite direction. This boustrophedonic motion is evident in the data 
records as alternating signs on the X spacing entry. 
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Table 1. Data File Structure 


Header Record, Once /Scan 

Header Record Code "H" A1 

Metric/English "M" or "E" A1 

Probe Calibration A6 

Module ID A2 

Test Number 14 

Run ID 12 

Scan ID Number 13 

Software Update Number 14 

Hour 12 

Minute 12 

Second 12 

Month 12 

Day 12 

Year 12 

Probe Type II 

Scan Type "1" = Rectilinear II 

Channel 0 Amp Ratio Code II 

Channel 1 Amp Ratio Code II 

Channel 2 Amp Ratio Code II 

Channel 3 Amp Ratio Code II 

Scale Factor Code II 

Number of Raster Repeats 12 

Spacing X F6.2 

Raster X Delta F6.2 

Raster Y Delta F6.2 

Initial X Position F6.2 

Initial Y Position F6.2 

Initial Z Position F6.2 

Zero X Position F6.2 

Zero Y Position F6.2 

Zero Z Position F6.2 
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Table 1. Data File Structure (Cont'd) 


Header Record, Once/Scan 

Zero Data Value 

F6.2 

Reference Intensity Value 

F6.2 

Local Intensity Value 

F7.2 


Data Record, Once/Line of Scan 


Data Record Code "D" Al 

Number of Data Points 13 

Hour 12 

Minute 12 

Second 12 

Spacing X and Direction (sign) F6.2 

Reference Intensity at Time HH:MM:SS F6.2 

X Position of Data Point 1 F6.2 

Y Position of Data Point 1 F6.2 

Data Point 1 F6.2 

Data Point 2, etc. F6.2 


End Record, Once /Scan 


End of Scan Code "E" 


Al 

Reference Intensity 


F6.2 

Check Sum of Absolute Value 

of 


All Data Points 


F12.2 

Summation Value (Negative = 

Overflow) 

19 

Points Summed 


14 

Number of Lines 


13 

End Hour 


12 

End Minute 


12 

End Second 


12 

End of Scan Character "!" 


Al 
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All entries in all the types of records were separated by a comma. All 
numerical entries were ASCII characters, with or without a leading positive or 
negative sign. All floating point numbers required a decimal point in the 
entry while integers did not. The only character that was not from the above 
list was the last character of the scan. This character was an ASCII 

code octal number 41. 

The acquisition software accessed the data port and verified the 
character type and structure of each record. If both were correct, then the 
information was decoded from ASCII to binary numbers when appropriate and 
stored on magnetic tape in a binary data file. This procedure was to be 
accomplished at the end of each complete raster of the flux mapper. At that 
time, the flux-mapper memory would be initialized and the next raster would be 
taken . 


The initial version of the data-output program was designed to read the 
binary data file on magnetic tape and print out the numbers as they were 
gathered from the flux mapper. No initial processing other than formatting of 
the data was to take place. 


D. PROBLEMS ENCOUNTERED 

Several immediate problems were encountered using these initial versions 
of the software, including polling with the operating system and transmission. 
The latter problem was never solved. 

As mentioned previously, the polling problem with the operating system 
placed a severe restriction on data acquisition. During the polling of all 
serial inputs, the operating system stores characters in a buffer for only the 
serial input being polled. The other inputs are ignored. The flux-mapper 
controller transmitted data to the serial port at a regular pattern regardless 
of the status lines in the RS-232C cable. This regular transmission pattern 
occasionally overlapped with the operating system's polling of the other 
ports. The result was that the character transmitted at that time was lost. 

To overcome this problem, the priority of the data acquisition task was 
altered . 

The RSX-llM operating system features a series of priority levels at 
which tasks can be assigned. All tasks with the same priority are polled in a 
round-robin fashion, and all of the users with the same priority level are 
given an equal opportunity to use the resources of the central processor unit 
(CPU) . Normally, tasks are given a priority level of 50 out of a maximum of 
250. Some tasks that require more of the CPU's resources are assigned higher 
priorities. An example of this is the text editor. It is normally 
installed at a priority of 65 because its interactive nature requires more of 
the CPU's resources. 

The data-acquis ition program required all of the CPU's resources for 
recording transmitted data from the flux-mapper controller. The acquisition 
program was installed at a high priority, 249, prior to execution. At 
execution, this task occupied virtually all of the CPU's resources and, 
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therefore, gathered all the data without difficulty. However, monopolizing 
all of the resources essentially rendered the multi-user system a single-user 
system. 


During initial use of the acquisition software, it was found that decode 
errors were occurring in the data string. These errors occurred randomly in 
the string, but in each occurrence the program would abnormally exit and 
abnormally close the binary data file on magnetic tape as a file of zero 
length. This was particularly annoying when the program would fail and all 
but the last few scans had been transmitted. In addition, this zero-length 
file on magnetic tape presented problems of playback of subsequently recorded 
data and also used a file name for null data. 

It was suggested that perhaps the flux-mapper controller had transmitted 
the incorrect character, causing the program to fail. This was checked by 
recording the output of the flux-mapper controller onto a digital data 
cassette tape unit. When several runs of the same set of data were recorded 
onto the cassette, a direct playback showed that these errors did not occur in 
every run nor at the same location of the run in which they occurred. 

It was decided that, to minimize the procedure for the personnel at the 
PDTS, the option of the analysis and printout immediately after a test was 
abandoned. Instead, the acquisition of the data was ensured in a simple way. 
The method used was the digital data cassette recorder mentioned above. At 
the end of each raster scan, the PDTS personnel transmitted the data from the 
flux-mapper controller to the digital data cassette three times. The cassette 
was then sent to JPL for archiving and printout. The three runs ensured that 
at least one complete and correct run had been stored. 

The digital data cassette recorder presented its own set of problems. 

In the record mode, the cassette fills a buffer and then transfers the data to 
tape in a single block. If the recorder was switched from the record mode, 
the data in a partially filled buffer was not transferred to tape; it was 
lost. Therefore, the recorder, once set to record mode, was left in record 
mode until all rasters of all the tests had been recorded for that particular 
day. The recorder presented a problem when tapes were recorded on 
one side and then reversed. Data on the first side appeared to be erased. 

The staff at the PDTS was instructed not to reverse Che tapes. With these 
instructions, the data was archived and printed at JPL at a later time. 

Finally, the use of the digital data cassette recorder allowed the 
PDP-11/34A to be used to acquire data from the Acurex Autodata-Nine data 
logger during the flux— mapper runs. The data-acquisition program monitored 
the data for various warnings, such as low cooling water flow, and displayed 
this information on a monitor. The data-acquisition program for the data 
loggers and for the flux-mapper controller are mutually exclusive tasks 
because of the nature of the operating system. It was felt that the warning 
alarms were an additional bonus. 
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SECTION III 


LOGGING METHOD USED 


A. DESIGN PHILOSOPHY 

The software was designed with three general design criteria in mind: 
compatibility with the flux-mapper controller, ease of use in the field, and 
error-checking capability. A large trade-off resulted for the criterion of 
ease-of-f ield use; thus, processing the data was made much more difficult. 

B. AVAILABLE HARDWARE 

The digital data cassette recorders, originally purchased as a backup 
unit for the data loggers, proved to be the important link between the 
flux-mapper controller and the PDP-11/34A. Two recorders, one at the test 
site for the recording of data and one at JPL for playback of data, were 
available. The one at the test site was set to the output characteristics of 
the flux mapper. The unit at JPL for data replay was set to an available 
port. The tapes used were cassette tapes of digital computer quality that 
contained archived as well as processed data. 


C. SOFTWARE WRITTEN 

The use of the digital data cassette recorders altered the original 
software tasks. The software tasks were subsequently divided into three 
parts: (1) software to read the data from cassette tape and write onto a disk 
storage medium, (2) software to read from the disk storage medium and transfer 
to a nine-track magnetic tape for archiving, and (3) software to print out the 
results from magnetic tape. 

The program to read the data from cassette tape and place it onto disk 
storage media was called CASTAP.FTN. This program had to be an installed task 
with logical unit number 3 reassigned to the input port. This task also 
required a high priority to bypass the polling option. (See the RSX-llM V.3.1 
operator's manual for details.) 

The program CASTAP transferred a given raster scan from the data 
cassette to a file on disk. The actual ASCII characters were transferred — 
no conversion of any type was made on the data elements. The program allowed 

the entry of the disk file name and allowed the selection of the scan to be 

stored on disk. The scan numbering system started with the present scan and 
incremented each time that the end of scan character, the " t", was found. At 

the conclusion of the data transfer to disk, the cassette recorder was turned 

off by the program. 
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The next step was to visually check the data file for incorrect 
characters. This was a rather poor method of error checking, but considering 
the limits of time and personnel available, it was the only one possible. The 
most common extraneous characters found in a record were lower-case characters 
and carriage returns embedded in data records. On finding errors, one of two 
options was available. The entire scan could be rerecorded, using one of the 
other three scans recorded; or a visual inspection of the data cassette 
rasters, using the data cassette recorder playing back directly into a 
terminal, would display the data from one of the other scans. The data in the 
disk file could then be changed to the correct value by using the text 
editor. To ensure that the file was examined, a comma, had to be added 

to the end of the disk data file following the end of raster character, the 
"!". If this comma were not added, then the file was deemed incomplete, and 
the archiving program would abort abnormally. 

The program to transfer the data from disk to a binary file on 
nine-track magnetic tape was called FMPCAS.FTN. This program read the ASCII 
file on disk, checked the record structure, decoded the ASCII to binary, and 
stored the binary in a nine-track magnetic tape to be initialized, if new, and 
be software— mounted . The magnetic tape file was opened as a Fortran logical 
unit and interfaced with the operating system. 

Once the program had been transferred to magnetic tape, the data could 
be printed out onto the line printer using the program FMPRINT.FTN. This 
program required the input of the data file name as well as requiring that the 
magnetic tape be mounted. Each file had to be called separately for printing; 
however, several copies could be produced with each call of this program. 


D. PROBLEMS ENCOUNTERED 

The problems with the software were quite evident. The delay between 
acquisition at the PDTS and the final printout, the slowness of the 
acquisition itself, and the tedious data scanning using the editor were all 
problems that would have been attacked, had there been time and personnel 
available. Twice a year, perhaps, rasters were taken over the course of a 
two-week period. This operation, therefore, was relatively infrequent and not 
commanding priority of time and personnel. 


E. SOFTWARE UPDATES 

From this version of the acquisition software, a major alteration to the 
data format was executed by the flux-mapper-controller programming group. 
Engineers analyzing the flux-mapper data requested time information for each 
scan. Because a typical complete raster would normally take from one-half 
hour to one-and-one-half hours, it was felt that these data would be very 
important in correlating the flux-mapper intensities to weather data such as 
insolation. The programmers added time information at the beginning of each 
data record. In addition, the output order was changed to its present 
boustrophedonic form. Both these changes occurred in June 1981. The software 
used for data acquisition included these changes. The plotting software (see 
below) included both formats. 
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SECTION IV 


DISPLAY SOFTWARE 


A. DESIGN PHILOSOPHY 

It was decided that a three-dimensional plot of the flux-mapper 
information would be the best way to represent the data for quick review. 

(See Appendix B.) It was felt that the plots should be easily understandable. 
Flexibility in plotting was also considered an important criterion. Viewing 
the flux map from various angles would greatly aid in understanding the result- 
ing patterns. The plots were designed to provide the data in a uniform manner 
to allow easy comparison between rasters of different distances from the focal 
plane. 


B. AVAILABLE HARDWARE 

The same basic minicomputer hardware was available. For printing and 
plotting, a Versatec Model 1100 electrostatic printer/plotter was used. The 
associated software to interact with the plotter portion of the device was 
purchased also from Versatec specifically for RSX-11M version 3.1. 


C. SOFTWARE WRITTEN 

The display software was a set of programs to input, reformat, 
calculate, and plot the flux-mapper raster data into a three-dimensional plot 
viewed from any location. The software was modeled after a plotting package 
in use on the Univac 1100 computer at the time of this development. Because 
of the size of the task involved, three -separate major programs were written, 
along with several other minor- programs . 

The first— minor program, TPDK.FTN, read the data from the nine-track 
magnetic -tape and transferred it to a general file on disk in the appropriate 
format. The same general file name on disk was used each time the program was 
called. Because' the program accessed the tape file as a Fortran logical unit, 
the magnetic tape had to be software-mounted. 

The first major program, GPL0T1D.FTN, read the data from the disk file 
and reformatted it for use by the other plotting routines. This routine also 
determined the type of plot desired, such as which view and the presence of 
contour lines. 

The second major program, GPL0T2D.FTN, performed the actual 
three-dimensional calculations, including the hidden line algorithm. The 
results were then placed in a file for final plotting in the last phase. 
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The third major program, PLOTERD.FTN, took the output from the previous 
file and created the Versatec plotter file. In this program, all of the 
interfacing to the Versatec software was included. A short, 30-character 
title was also requested as well as changes in scale factor. The resulting 
files, VECTRl.BIN and PARM.BIN, contain the plotting instructions for the 
Versatec plotter. 

The last step was to invoke the plotter, using the Versatec-supplied 
program, RASM.TSK. This program accessed VECTRl.BIN and PARM.BIN to produce 
the actual plots. 

A set of modified versions of the major plotting programs was created 
for use at the PDTS. The number of optional features, such as viewing 
location, were fixed. These were denoted as FPL0T1.FTN, FPL0T2.FTN, and 
PLOTER.FTN. 

A second set of modified versions was created to view the plot from 
directly overhead, namely, a contour plot showing isointensity lines of solar 
radiation. The same general steps were followed in this set of programs, 
named CPLOT.FTN, CDRIVE.FTN, and SUMRAD.FTN. 

These plots, combined with the data printouts, provided a complete 
picture of the solar irradiation. 
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SECTION V 


RECOMMENDATIONS AND CONCLUSIONS 


A. SOFTWARE RECOMMENDATIONS 

The software recommendations discussed here would have been Implemented, 
given sufficient time and personnel. Because of the software's infrequent 
use, most of the changes and modifications were made either just prior to or 
immediately after its use, when its priority was high. 

The plots and data printouts were analyzed and used by JPL personnel. 

The filing system that was used did not match the one used by the data- 
acquisition software. Much of the correlation of older data has been from 
descriptive information included on the plots and printouts. It would have 
been very useful to provide more space for descriptive information on both the 
plots and printouts. 

The operating system and the problem with the system polling may have 
required some system programming modification. At the time, no available 
personnel had the expertise to examine the problem. Also, the operating 
system was an older, unsupported version. The newer versions may, in fact, 
alleviate or mitigate this buffer/polling problem. 

The tedius work of going through the disk file after transfer from the 
digital data cassette using the text editor might have been alleviated by a 
program that would scan the ASCII data for inappropriate characters. Time did 
not permit the writing of this program. 

Finally, with the correct type of hardware, it should be possible to 
reproduce the three-dimensional flux-mapper plots on a video terminal equipped 
for graphics output. With interactive features, this would allow for easier 
interpretation of data plots. 


B. HARDWARE RECOMMENDATIONS 

During the initial checkout of the software, the unavailability of the 
flux-mapper hardware made correction and modification difficult. It would 
have been helpful if a hardware simulator had been available. An alternative 
would have been to perform all of the development work, both hardware and 
software, in one location. 

The Versatec printer/plotter was an older model, no longer in production. 
It normally required several minutes to produce a plot. A faster plotter 
would have been useful. 
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c. 


CONCLUSIONS 


The flux-mapper software, as the flux mapper itself, was a 
tool. Although the use of the software was cumbersome at times, 
consistently provided useful plots and printouts to evaluate the 
of point-focusing parabolic dishes . 


laboratory 

it 

performance 
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APPENDIX A 


LIST OF COMPUTER PROGRAMS 


This Appendix contains one of the representative sets of programs used 
to analyze and plot the flux-mapper data. The set included herein is the 
standard set used for analysis of most flux-mapper runs. This software was 
available to the PDTS staff; analysis of data contained on digital cassettes 
was usually performed at the main JPL facility. 


The programs included in this representative set of software are 


FMTPDK.FTN 

FPLOT1.FTN 

SET32.FTN 

CLSET.FTN 

FPL0T2.FTN 

SEG.FTN 

CTCELL.FTN 

DRAW . FTN 

PLOTER.FTN 


Flux-mapper mag tape to disk program 

First stage of the field three-dimensional plot system 
Subroutine for three-space to two-space transformation 
Subroutine to calculate contour values 

Second stage of the field three-dimensional plot system 
Subroutine to write plotting data into a file 
Subroutine to compute contour lines 

Subroutine to draw visible part of a line between two 
points 

Third stage of the field three-dimensional plot system 


All of the software generated for the flux-mapper task is not included 
here due to space constraints. Copies of this software were transferred to 
Sandia National Laboratories-Albuquerque (SNLA) and are available through that 
organization. 


Major software packages developed for the flux mapper but not listed in 
this Appendix include the following: 


CPLOT.CMD 
CPLOT.FTN 
CPLOTO. FTN 

CDRIVE.FTN 
GPLOT . CMD 
GPLOTlD. FTN 
GPL0T2D. FTN 

PLOTERD . FTN 
SUMEVN.FTN 
SUMP RT. FTN 
SUMRAD. FTN 
SUMRADO. FTN 

SUMUP.FTN 


Command file to run contour plot software system 
First stage of contour plot system 

First stage of contour plot system for old flux-mapper 
files 

Second stage of contour plot system 

Command file to run generalized plot software system 
First stage of generalized three-dimensional plot system 
Second stage of generalized three-dimensional plot 
sys tern 

Third stage of generalized three-dimensional plot system 
Superimposes two flux-mapper files of unequal size 
Prints flux-mapper data to printer 
Integrates flux over a user-defined area 
Integrates flux over a user-defined area for old 
flux-mapper files 

Adds the contents of two flux-maper files together 


All of the above software was available for use by the PDTS staff. However, 
only the standard field packages included in this Appendix were consistently 
used . 
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C FMTPDK.001 

C 

C Flux mapper tape to disk program. 

C This program is designed to 

C dump the contents of the flux mapper 

C tape onto disk for plotting 

C and other uses. 

C 

C JPL PFDRT. Written by Stephen Pitchie 

C 

C .001 15-SEP-80 INITIAL VERSION 

C 

C LINK: Uses standard libraries 

C 


£•**** ******************************************* 

LOGICAL* 1 IY, YY, ICD, FILE { 14 ) , BK{ 14 ) , IN , 

* NE 

DIMENSION DTA(64),ISCCM(41) < FDTA(M) 

INTEGER 1M( 3 ) , Z, Y, HER (9 ) , STYPE, DAY, YEAR, RID, TPRQRE, PCAL{ 3 ) , SID, 
*END, SUNO, CSF ( 4 ) , SFC , TNO 


C 


20 


23 


REAL LIV 
DATA YY/'Y'/ 

DATA FILE/ 'M' , 'T', '0', * : ' , 'F' , 'M* , ’0' , 'O', *0', '0* 
**A* , 'T'/ 

DATA IHR/ 'FH'/» IDR/ ' FD ' / , ITR/ 'YE' / 

ITAPE = 4 


' , 'D' , 


TYPE Enter Flux mapper file name in XXXX.YYY format' 

ACCEPT 902, (FILE(I ) , 1=5, 14) 

TYPE 903, FILE 

CALL ASSIGN (ITAPE, FILL, 14) 

REWIND ITAPE 
NO = 1 

OPEN (UNIT=1,NAME=’ TEMP. DAP’ ,TYPE='NEW' , FORM= 1 UNFORMATTED 1 , 
*ACCESS= ' DIRECT ' , ASSOCIATEVARIABLE=NO, RECORDS IZE=100, DISP= ' SAVE ' ) 
READ (ITAPE) 

READ (ITAPE) 

READ (ITAPE) 

READ (ITAPE) 

READ (ITAPE) (ISCCM(I), 1=1, 10) 

WRITE (1 'NO) (ISCCM(I), 1=1,10) 

READ ( ITAPE )IH 
IF (IH.EQ. IHR) GO TO 20 


TYPE 905,111 
GO TO 50 

READ ( ITAPE ) METRIC , PCAL, MODID , TNO, RID, SID, SUNO , TM , MCNTH , DAY, 
*YEAR, TPROBE, STYPE , CSF, SFC, NRR, SX, DM, DY, XI , YT , ZI , ZX, ZY, ZZ , ZDV, RIV, 
*LIV 

WRITE ( 1 * NO ) IH , METRIC , PCAL, MODID , TNC , RID , SID , SUNO , TM , MONTH , DAY, 
♦YEAR, TPROBE, STYPE, CSF, SFC, NRR, SX, DX, DY, XI , YI , ZI , ZX, ZY, ZZ, ZDV, RIV, 
*LIV 

READ (ITAPE) ID 
IF (ID.EQ.IDR) GO TO 25 
IF (ID.EQ.ITR) GO TO 29 
TYPE 910, ID 
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GO TO 50 

25 READ { ITAPE ) N, II, 12, 13, (DTA(I) , 1=1, N) 

C IF (DTA(l).GE.O.O) GO TO 27 

GO TO 27 

C RE-ORDERING INPUT DATA 

DO 26 1=5, N 
JJ =65-1 
RDTA(JJ) = CfTA(I) 

26 CONTINUE 

WRITE(1 ’NO) ID,N,dta(l),dta(2),dta(3),dta(4), (RDTA(l ) , I=JJ, 60) 
GO TO 20 

27 WRITE (1 ’NO) ID,N, (DTA(I ) , 1=1 ,N) 

2P CONTINUE 

GO TO 23 
29 CONTINUE 

READ ( ITAPE ) RIN , DCS , SV , NPS , NLINES , IER, IEM,IES 
WRITE (1 'NO) ID.RIN, SCS,SV, NPS, NLINES, IER, IER, IER 
CLOSE (UNIT=1) 

CALL CLOSE (ITAPE) 

50 STOP 

901 FORMAT (Al) 

902 FORMAT (20A1) 

903 FORMAT (111 10010 ) 

904 FORMAT (1H F7.2) 

905 FORMAT (' HEADER RECORD DOES NOT MATCH’, A2) 

910 FORI AT (' DATA RECORD DOES NOT MATCH' ,A2) 

END 



PROGRAM FFLCT1 

LOGICAL* 1 YY,DILE(10),CDILE(8),IY,BELL 
DIMfNSION X(33 ) ,Y(33) ,Z(1089) ,M(2178) , 

S(6) ,DTA(40) ,Y1(33) 

DIMENSION MXS (2 ) ,MXF(2) ,MXJ(2) ,MYS(2) ,MYF(2) ,MYJ( 



CCMCN /SRFBLK/ LIMU(102.4) ,LIML(1024) 

,CL(41) ,NCL , 

* 

LL , FACT , IROT 

,NDRZ 

★ 

NUPPER , NRSWT , BIGD 

, EMIN 

* 

UMAX , VMIN , VMAX 

, RZERO , 

★ 

NOFFP , NSPVAL ,SPV 

, BIGEST 


CCMMCN /PWRZ1/ XXMIN ,XXMAX , YYMIN , YYMAX 

★ 

ZZMIN , ZZMAX , DELCRT 

, EYEX , 

it 

EYEY , EYEZ 



CCM-1CN /SET/CVAL(10) , NCHAR , IHEAD(IO) ,XDIST,NCN ,NVAL 
CCJ-WCN /SETN/ NN 
DATA BELL/'*7/ 

DATA SFVAL,IOFFP/0.0,0/,YY/'Y'/ 

DATA STEREO /O.O/ 

DATA IFR, ISTP , IRCTS, IDRX, IDRY, IDRZ , TOPPER, ISKIRT, NCLA/ 

* 1 , 0 , 0 , 1 , 1 , 0 , 0 , 0 , 6 / 

DATA THETA, HSKIRT, CHI , CLO,CINC/ 

* . 02 , 0 ., 0 ., 0 ., 0 ./ 

DATA S/40.,10.,10. ,0.,0.,0. /MDZ, NX, NY,NCHAR/33, 33, 33, 10/ 

NRSWT = 0 
IEND = -9999 
ISPVAL = -999 
BIGEST = 1.F.37 
END = -9999. 

NO = 1 

TYPE *, ' THIS PROGRAM PERFORMS A 3-D PLOT OF FLUX MAPPER DATA* 
OPEN ( UNIT=4, EAME=‘ TEMP. CAP ' , TYPE=’ OLD ' ,FORM=’ UNFORMATTED ' , 
1ACCESS=' DIRECT' ,ASS0CIATE\M13AELF,^0, RECORDS IZE=100) 

905 FORMAT (10A1) 

906 FORMAT (RAl) 

30 CONTINUE 

TYPE *, ' Enter nimber of X points in scan, INTEGER' 

ACCEPT 911, NX 

TYPE Enter number of Y lines in scan, INTEGER' 

ACCEPT 911, NY 
911 FORMAT (15) 

907 FORMAT (10A2) 

TYPE Do you want Front(l) or F,ide(2) view? Biter 1 or 
* 2 INTEGER' 

ACCEPT 911, NV 
IF (NV.EQ.l) GO TO 304 
S(l) = 10. 

S(2) = 40. 

304 CONTINUE 

908 FORMAT (F7.2) 

305 TYPE *, 'Do you want contours? Y(es) or N(o)’ 

ACCEPT 905, IY 

IF (IY.GT.8288) IY = IY - 32 
IF (IY.EQ.YY) IDRZ = 1 
READ (4 'NO) 

READ (4'N0)IH,MM,M1,M2,M3,M4,M5,M6,M7,M8,M9,N,N1,N2,N3,N4,N5,N6, 
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* N7,NG,N9, 11, 12, I3,EDX,X1,X2,X(1),X3,X4,X5,X6,X7,X8,X9,X10 
IF (I2.EQ.0) AJJ = .0125 
IF (I2.E0.1) AIJ = ,5 
IF (12. EC-?) AIJ = 1. 

IF (12. EC- 3) AIJ = 2. 

IF (I2.EQ.4 ) AIJ = 5. 

IF (I2.EQ.5) AIJ = 10. 

IF (I2.EQ.6) AIJ = 20. 

IF (12. EC. 7) AIJ = 50. 

IF (I2.EQ.F) AIJ = 100. 

IF (I2.E0.9) AIJ = 200. 

CO 31 1=2, MX 
X(I) = X(I-l) + DCX 

31 CONTINUE 
X1G = X(l) 

X2G = X(NX) 

KSWT = 0 
KFLG = 0 

TYPE *,' IS DATA PRE-JUNE 1901? (y/N)‘ ' 

ACCEPT 905, IY 
IF (IY.GT.9G) IY = IY -32 
IF (IY.NE.YY) KFLG = 1 
DO 3G K=1,NY 

READ (4 'NO)IC,N, (DTA(l) , I=1,N) 

Y1(K) = DTA(4) 

IF (N.EQ.NX+4) GO TO 32 
KCK = N 
N = NX + 4 

TYPE *,' ***** WARNING, NOT ENOUGH X VALUES ****** 

JKX = KCK - 4 
TYPE 917, JKX 

917 FORMAT (* ONLY*, 13, * VALUES FOR X’) 

TYPE 91 R, BELL 
TYPE 918, BELL 
91 P FORI'IAT (1H Al) 

32 CONTINUE 

IF (KFLG. EG • 1 ) GO TO 34 
DO 33 J=5, N 

Z( J-4-*NX*(K-l ) ) = DTA(J) 

ZCMAX = AMAXl ( ZCMAX, DTA( J ) ) 

33 CONTINUE 
GO TO 38 

34 IF (KSWT.NE.O) GO TO 36 
DO 35 J=5,N 

Z ( J -4+NX* (K-l ) ) = DTA(J) 

ZCMAX - AMAXl ( ZCMAX, ETA( J ) ) 

35 CONTINUE 
KSWT = 1 
GO TO 38 

36 DO 37 J=5,N 
Z(J-4+NX*(K-l)) = DTA(N+5-J) 

ZCMAX = AMAXl (ZCMAX, DTA(N+5-J ) ) 

37 CONTINUE 
KSWT = 0 

38 CONTINUE 
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IF (KFLG .EQ. ] ) GO TO 40 

EC 39 J=1,NY 

Y(J ) = Y1 (NY^T+1 ) 

39 CONTINUE 
GO TO 42 

40 DO 41 J=1,NY 
Y(J) = YJ (J) 

41 CaJTINUE 

4? CONTINUE 

CLOSE (UNIT=4) 

C PRINT 9C1,X,Y,Z 

Y1G = Y(l) 

Y2G = Y(NY) 
za-iAX = zcmax*aij 

901 FORMAT (1H 9EI4.7) 

NON = 1 
MMXX - MDZ 
NNXX » NX 
NNYY = NY 
STER = STEREO 
NXPl = NNXX + 1 
NYPl = NNYY + 1 
NLA = NCIA 
NSPVAL = ISPVAL 
NOFFP = IOFFP 
SFV = RPVAL 
NDR'Z = IDRZ 
IF (IDRZ.NE.O) 

* CALL CLELT ( Z , MMXX , WXX , MIYY , CHI , CIO , CINC , NLA ,40, CL, NCL , 

* ICNST, NOFFP, SPV, DIGEST) 

IF (IDFZ.NE.P) NDRZ = 1 - ICNST 
STHETA = SIN(STER*THETA) 

CTHETA = COS ( RTER*THETA ) 

RX = S(l) - S (4 ) 

RY = S (2 ) - S (5 ) 

RZ = S{3 ) - S (6 ) 

D1 = SORT ( RX*RX+RY* RY+RZ* RZ ) 

D2 = SORT ( RX* RX+RY* RY ) 

DX = 0. 

DY = 0. 

IF ( STEREO. EQ. 0. ) GO TO 102 
D1 = Dl*STEFEO*THETA 
IF (D2.OT.0.) GO TO 101 
DX = Dl 
GO TO 102 

101 AGL = ATAN{FX/-RY) 

DX = Dl TOPS (AGL) 

DY ~ Dl*SIN (AGL) 

102 IROT = IROTS 
NPIC = 1 

IF (STER.NE. 0. ) NPIC = 2 
FACT = 1. 

IF (NRSWT.NE.O) FACT = RZERO/D1 

IF ( TOTP.EO. 0. AND. STER.NE. 0. ) IROT = 1 

XDIST = .5*( (1024. /102. )-. ]P*FLOAT(NCHAR) ) 
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IP1C = 1 

NUPPER = I UPPER 
SIGN1 = I PIC*? - 3 
EYEX = S(l) + SIGN1*DX 
POIX = S(4) + SIGN1*DX 
EYEY = S(2) + SIGNl*DX 
POIY = S(5) + STGN1*DX 
EYEZ = S(3) 

POIZ = s(r>) 

LL = 0 

GYLL SET32 (POIX, POIY, POIZ, EYEX, EYEY, EYEZ, 3 ) 

LL = IPIC + 2*ISTP + 3 
IF (STER.EQ.O.) IX = 1 
IF (NRSWT.NE.O) GO TO 107 
XXMIM = X(l) 

XXMAX - X(NNXX) 

YYHIN = Y(l) 

YYMAX = Y(NNYY) 

UMIN = BIGEST 
VMIN = BIGEST 
ZZMIN = BIGEST 
UMAX = - UMIN 
VMAX = -VMIN 
ZZMAX = — ZZMIN 

C PRINT 901, UMIN, UMAX, VMIN, VMAX, ZZMIN, ZZMAX 

DO 104 J=1,HJVY 

DO 103 I=1,NNXX 

ZZ = Z( I+NX* ( J-l ) ) 

IF (NOFFP.EQ. 1 .AND.ZZ.Fp.SPV) GOTO 103 
ZZMAX = AMAX1( ZZMAX, ZZ) 

ZZMIN - AMIN1( ZZMIN, ZZ) 

CALL SET32 (X< I ) ,Y(J) ,Z(I-fNX*( J-l)) ,UT, VT,DU, 2) 
UMAX = A MAXI (UMAX, UT ) 

UMIN = AM INI (UMIN, UT) 

VMAX = AMAX1 (VMAX.VT) 

VMIN = AMIN1 (VMIN, VT) 

103 CONTINUE 

104 CONTINUE 

C PRINT 901, UMIN, IMAX, VMIN, VMAX, ZZMIN, ZZMAX 

WIDTH = UMAX -UMIN 
RIGHT = VMAX -VMIN 
DIF = .5*(WJDTHHHIGHT) 

C PRINT 901, WIDTH, I1IGHT, DIF 

IF (DIF) 105,107,106 

105 UMIN = UMIN + DIF 

UMAX = UMAX - DIF 
GO TO 107 

106 VMIN = VMIN - DIF 

VMAX = VMAX + DIF 

107 GALL SET32 (POIX, POIY, POIZ, EYEX, EYEY, EYEZ, 1) 

C TYPE 900,mXX,NNYY 

900 FORMAT (1H 13010) 

DO 109 J=1,NNYY 

DO 100 I=1,NNXX 

CALL SET32 (X( 1 ) ,Y( J ) , Z( I+NNXX* ( J-l ) ) ,UT, VT, DU, 2) 
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910 

103 

109 

110 

111 

112 

113 

114 

115 

116 

117 

113 

119 


M( 2*1-1 +2*NNXX*(J-1)) = UT 
M(2*I+2*NNXX*(J-1) ) = VT 
PRINT 910, I, J,M(1, I, J) ,1, J,M(2, I, J) 

FORMAT (III 2(213, 2X, 15, 7X)) 

PRINT 900, M(l, I,J),M(2,1,J) 

PRINT 901,UT,VT 

CONTINUE 
CONTINUE 
DO 110 K=1 , 1 024 
LIMU(K) = 0 
LIML(K) = 1024 
CONTINUE 
NXPASS = 1 

IF (S (1 ) .CE.X(NNXX) ) GO TO 113 
IF (S{1) .LE.X(l)) GO TO 114 
DO 111 1=2, NNXX 
LX = I 

IF (S(l ) .LE.X(I)) GO TO 112 
CONTINUE 
MXS(l) = LX-1 
MXJ(l) = -1 
MXF{1) = 1 
MXS(2) = LX 
MXJ(2) = 1 
MXF(2) = NNXX 
NXPASS = 2 
GO TO 115 
MXS (1 ) = NNXX 
MXJ(l) = -1 
MXF(l) = 1 
GO TO 115 
MXS(l) = 1 
MXJ(l) = 1 
MXF{ 1 ) = NNXX 
NYPASS - 1 

IF { S ( 2 ) • GE . Y { NNYY ) ) GO TO 113 
IF (S(2) .LE.Y(l) ) GO TO 119 
DO Ilf J=2 , NNYY 
LY = J 

IF (S(2).tE.(Y(J)+.02)) GO TO 117 
CONTINUE 
MYS(l) = LY - 1 
MYJ(l) = -1 
MYF(l) = 1 
t IYS ( 2 ) = LY 
MYJ (2) = 1 
MYF(2) = NNYY 
NYPASS = 2 
GO TO 120 
MYS(l) = NNYY 
MYJ{1) = -1 
MYF(l) = 1 
GO TO 120 
MYS(l) = 1 
MYJ ( 1 ) = 1 
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MYF(l) - NNYY 

120 NN = 1 
C PRINT 901, X,Y 

C PRINT 901, Z 

C PRINT 901 , XXMIN, XXMAX, YYMIN, YYMAX, ZZMIN, ZZMAX, PELCRT 

C PRINT 901, EYEX,EYEY,EYEZ, FACT, EIGD, DIGEST, RZERO 

C PRINT 901, EMIN, UMAX, VMIN,VMAX 

C PRINT 900, HIEAD, NCL, II, , I ROT, NDRZ , NUPPER, NRSWT, NOFFP, NSPVAL 

OPEN ( UNIT=3 , NAME= ’ SCRT . DAT ' ,TYPE='NEW' , FORM= ' UNFORMATTED ' , 

* ACCESS= ' DIRECT 1 ,ASSOCIATLVARIABLE=NN, RECORDS IZE=1 090) 

WRITE (3 'NN )IHEAD,X,Y, XXMIN, XXMAX, YYMIN, YYMAX, ZZMIN, ZZMAX, 

*PELCRT, EYEX, EYEY, EYEZ , NCL, LL,FACT, TROT , NDRZ , NUPPER , NRSWT , 
*DIGD,UMIN, UMAX, VMIN, VMAX, RZERD, NOFFP, NSPVAL, £PV, DIGEST, CL 

* , NX, NY, MDZ , X1G , X2G , Y3 G , Y2G , ZCMAX, AI J 

C PRINT 900, LIMU 

WRITE ( 3 ' NN ) ( LIMU ( I ) , 1=1 , 1024 ) 

C PRINT 900, LIML 

VJRITE (3'NM)(LiriL(l), 1=1, 1024) 

V TRITE (3'NN)(M(l), 1=1, 2178) 

WRITE (3'NN)(Z(I), 1=1,1089) 

PRINT 901 , PX, RY, HSKIRT, XDIST, CVAL, S, CL 

PRINT’ 900 , NXPASS , NYPASS , MXS , MXF , MXJ , MYS , MYF , MYJ , IDAY, NCHAR 
PRINT 900, NON, NW, 

WRITE ( 3 1 NN ) NXPASS , NYPASS , MXS , MXF , MXJ , MYS , MYF , MYJ , RX, RY, IDAY, 

HSKIRT, NDRZ , XDIST, CVAL, NCHAR, NON, NVAL, S, NNXX, NNYY, MMXX, ISKI 


CALL SET3? (X (1 ) , Y(1 ) ,Z(1 ) ,UT, VT.DUM, 3) 
CLOSE (UNIT=3 ) 

STOP 

END 
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SUBROUTINE SET32 (X, Y, Z , XT , YT , ZT, KFLAG ) 

C 

C Hi is routine implements the 3-space to 2 -space transformation 

C by Kuber, Szabo, and Giuliori, The Perspective Representation 

C of Rmctions of TUo VAriables, J. /CM 15, 2 193-204,1968 

C 

C X,Y,Z Are the 3-space coordinates of the intersection of the 

C line of sight and the image plane. This point can be 

C thought of as the point looked at. 

C XT, YT, ZT Axe the 3-space coordinates of the eye position. 

C 

C KFLAG — 2 arguments 
C 

C X, Y, Z Are the 3-space coordinated of a point to be 

C transfonned . 

C XT, YT 'flie results of the 3-space to 2-space transformation. 

C ZT Not used . 

C 

C If LL(in CCMMCN)=0, XT and YT are in the same scale and X,Y Z. 

C 

C The variable KFLAG has tvo possible variables 

C 1-oompute intersection of line of sight 

C 2-transform from 3-space to 2-spac l* 

C 

C NOTE.! 1 1 ! ! ! ! !J ! i ! 

C The KFLAG=3, 4 are special debugging flags and are not 

C part of the plot package. 

c 
C 

CCMMQN /FWFZ1/ XXMIN ,XXMAX , YYMIN , YYMAX 


* 

ZZMIN 

, ZZMAX 

, CELCRT , EYEX 

★ 

EYEY 

, EYEZ 



COMMON /SRFELK/ LIMU(1024) 

,LIML{1024 ) ,CL{41 ) ,NCL , 

★ 

LL 

,FACT 

, IROT , NDRZ 

★ 

NUPPER 

, NRSWT , BIGD ,UMIN 

* 

UMAX 

, VMIN 

,VMAX , RZERO 

★ 

NOFFP 

, NSPVAL ,SPV , BIGEST 


CCMMQN /SETN/ NN 




DIMENSION NLU(7) 

,NRU{7 ) 

, NB V ( 7 ) ,NTV(7) 


c 

C picture comer coordinates for LL=1 

C 

DATA NLU(1 ) , NRU ( 1 ) ,NBV( 1 ) ,NTV( 1 )/10, 1014, 10, 1014/ 
picture comer coordinates for IL=2 
DATA NLU(2),NRU(2),NBV(2),OTV(2)/10, 924,50, 964/ 
picture comer coordinates for LL=3 
DATA NLU(3 ) ,NRU(3 ) ,NBV(3 ) ,NTV(3 )/100, 1014,50, 964/ 
picture comer coordinates for LL=4 
DATA NLU(4) ,NRU(4),NBV(4) ,NTV(4)/10, 1014, 10, 1014/ 
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picture comer coordinates for LL=5 

DATA NLU{5 ) ,NRU(5) ,NBV(5) ,NIV{5)/10, 1014, 10, 1014/ 

picture corner coordinates for LL=6 

DATA NLU{f>) ,NRU(6) ,NBV(6) ,WTV{6)/10, 512,250,758/ 

picture comer coordinates for LL=7 

DATA NLU(7) ,NRU{7 ) ,NBV(7 ) ,NTV(7 )/512, 1014,256,750/ 

GO TO (1,2, 3, 4) KFLAG 
JUMP3 = 104 

IF (NOFFP.E0-1) JUMP3 = 103 

AX = X 

AY = Y 

AZ = Z 

EX = XT 

EY = YT 

EZ = ZT 

As much computation as possible is done during execution of 
SET32 since the transformation is called many times for each 
call to SFT32. 

DX = AX - EX 
DY = AY - EY 
DZ = AZ - EZ 

A more careful computation of direction Gosines. 

D = 0. 

T = AMAX1 (ABS(DX) ,ABS(DY) ,ABS(DZ) ) 

PRINT 901,AX,AY,AZ,EX,EY,EZ,DX,DY,DZ,T 
PRINT 900 , N0FFP 
FORMAT (1H 120010) 

FORMAT (1H 9E14.7) 

IF (T.EQ.0.0) GO TO 30 
PI = DX/T 

R2 = DY/T 
R3 = DZ/T 

D = S0RT(R1*R1 + R2*R2 + R3*R3) 

If D isn't ZERO 

COSAL = Rl/D 
COGBE = R2/D 
COSGA = R3/D 
D = D*T 
GO TO 40 

If D is ZERO, ray has no direction: assign direction down 
X-axis . 
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c 

30 CONTINUE 

COSAL = 1. 

COSUE =0. 

COSGA = 0. 

40 CONTINUE 

AL = ACCS (COSAL) 

BE = AC03 { OOSP.E ) 

GA = ACOS(OOSGA) 

SINGA = SIN(GA) 

C PRINT 901 , R1 , R2 , R3 , D, COSAL, COSBE, COSGA , AL, BE, GA, SINGA 

C PRINT 900, IL 

JUMP2 = 110 
IF (LL.EQ.O) GO TO 101 
JUMP2= 100 

DELCRT = NRU ( LL ) -NLU ( LL ) 

UO = UMIN 
VO = VMIN 
U1 = NLU(LL) 

VI = NBV(LL) 

U2 = NRU(LL) - NLU(LL) 

V2 - NTV(LL) - NBV(LL) 

IF(UMAX-UMIN) 52,51,52 

51 U3 = 0. 

GO TO 53 

52 U3 - U2/ ( UMAX-UMIN ) 

53 CONTINUE 

IF (VMAX-VMIN) 55,54,55 

54 V3 = 0. 

GO TO 56 

55 V3 = V2/ ( VMAX-VMIN ) 

56 CONTINUE 

U4 = NRU(LL) 

V4 = NTV(LL) 

C PRINT 901,TO,VO,U1,V1,U2,V2,U3,V3,U4,V4,UMAX,VMAX 

C PRINT 900, NRSWT 

IF (NRSWT. EC, 0) GO TO 101 
UO = -BIGD 
VO = -BIGD 
U3 = U2/(2.*BIGD) 

V3 = V2/(2.*BIGD) 

C 

C The 3-space point looked at is transformed inot (0,0) of the 

C 2-space. IVie 3-space Z axis is transformed into the 2-space Y 

C axis. If the line of sight is close to parallel to the 3-space 

C Z axis, the 3-space Y axis is chosen (instead of the 3-space 

C Z axis) to be transformed into the 2-space Y axix. 

C 

101 IF (SINGA.LT. 0.0001) GO TO 102 
R = 1./ SINGA 

JUMP = 105 

C PRINT 900, JUMP, JUMP2,JUMP3 

RETURN 

102 SINBE = SIN (EE) 

R = 1. /SINBE 
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c 

c 

c 

c 

2 


c 

c 


103 

104 


1111 

50 

C 

C 


105 

106 

107 

108 

109 

110 


C 

c 


3 


C 

C 

C 


4 


C 

C 

C 


JUMP = 106 

PRINT 900, JUMP, JUMP3 , JUMP2 
RETURN 


Transformation entry point 

XX = X 
YY = Y 
ZZ = 2 

PRINT 901,XX,YY, ZZ 
PRINT 900, JUMP3 
IF (JUMP3.EQ.104) GO TO 104 
IF (Z2.EQ.SPV) GO TO 109 

DENCM = (XX -EX ) * COSAL + (YY-EY)* COSBE + ( ZZ-EZ ) *COSGA 

if (dencm.ne.o.o) go to 1111 
Q = 1. 

CO TO 50 
0 = D/DENCM 
CONTINUE 

PRINT 901, DENCM, 0 

PRINT 900, JUMP 

IF ( JUMP. FC. 106) GO TO 106 

XX = ( (EX4Q*(XX-EX)-AX)*C05BE-(EY4Q*(YY-EY)-AY)*CCSAL)*P. 
YY = (EZ4Q*( ZZ-EZ )-AZ)*R 
GO TO 107 

XX = ( { EZ+Q* ( ZZ -EZ ) -AZ ) *G0SAL- ( EX+Q* { XX-EX ) -AX ) * COSGA ) *R 
YY = (EY+Q*(YY-EY)-AY)*R 
IF (JUMP2.EQ.110) GO TO 110 
XX = AMIN1(U4,AMA\1(U1,U14U3*(FACT*XX-U0))) 

YY = AM INI (V4 , AMAX1 (VI , V1+V3* { FACT*YY-VO) ) ) 

GO TO 110 
XX = NSPVAL 
YY = NSPVAL 
GO TO 110 
xr = xx 
YT = YY 

PRINT 901,XT,YT 
PRINT 900, JUMP2 
RETURN 
CONTINUE 

WRITE (3 ’NN)JUMP, JUMP2, JUMP3, EX, EY, EZ , COSAL, COSBE, COSGA, D, 
AX, AY, AZ, R,U0, U1,U2, U3,U4,VO,Vl, V2, V3, V4 
PRINT 900, JUMP, JUMP?, JUMP3 

PR IN'! 901 , EX, EY, EZ, COSAL, COSBE, COSGA, D, AX, AY, AZ, R 

PRINT 901, U0,U1 , U2, U3, U4, VO, VI, V2, V3, V4 

REIURN 

continue 

READ (4’NN) JUMP, JUMP2,JUMP3, EX, EY,EZ, COSAL, COSBE, COSGA, D, 
AX, AY, AZ , R, UO, U1 , U2 , U3 , U4, VO, VI , V2 , V3 , V4 
PRINT 900, JUMP, JUMP2,JUMP3 

PRINT 901 , EX, EY, EZ , COSAL, COSBE, COSGA, D, AX, AY, AZ , R 

PRINT 901, UO,Ul,U2,U3,U4,VO,Vl,V2,V3,V4 

RETURN 

END 
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FUNCTION ACOS (X) 

ACCS = l./£X2RT(l.+ATfiN(X)*ATAK(X)) 
END 


I 

\ 


f 

t 
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SUBROUTINE CLSET ( Z , MDZ , NO C, NY, CHI , CIO, CINC, NLA , NLM , CL, NCL, 
*ICNST,NOFFP,SPV, DIGEST) 

CCMMCN /SE1’/CVRL(10) ,MCHAR, IHFAD(IO) ,XDIST, KON,NVAL 
DIMENSION Z(10P9) ,CL(41) 

L0GICAL*1 GOODZ 
C 

C CLSET COMPUTES THE VALUES OF THE CONTOUR LEVELS 

C AND PUTS TRIM IN 

C CL 

C 

GOODZ = .FALSE. 

1CNST = 0 
GLO =CLO 
HA = CHI 
FANC = CINC 
CRAT = FLOAT (NLA) 

NCL = 0 

IF (HA -GLO ) 110,120,130 
110 GLO = HA 
RA = CID 
GO TO 130 
C 

C If HA and GLO are not set, set them to the values of DIGEST 

C and -DIGEST respectively 

C 

120 GLO = -BIC.EST 
RA = -GLO 
GO TO 140 

130 IF (FANC.EQ.O. ) FANC = (ILA-GLO)/(CPAT-l ) 

140 CONTINUE 

DO 150 J=1 , NY 

DO 145 1=1, NX 

IF (NOFFF. EC 1 - 1 . AND.Z( I+NX* ( J— 1 ) ) -EQ. SFV) GO TC 145 
IF (GOODZ) GO TO 146 
ZMIN - Z ( I+NX* ( J-l ) ) 

ZMAX = Z(I+NX*(J-I)) 

GOODZ = .TRUE. 

GO TO 145 

146 ZMIN = AMIN3 (Z( I+NX* (J-l)), ZMIN) 

ZMAX = AMAX1(Z( I+NX* (J-l )) ,ZMAX) 

145 CONTINUE 

150 CONTINUE 

Check HA and GLO to make sure they fall within the range of Z 
values being plotted; if not, set IIA=ZMAX and/or GLO=ZMIN. 
c 

IF (HA. GT. ZMAX) HA = ZMAX 
IF (GLO. LT. ZMIN) GLO = ZMIN 
C PRINT 901, HA, GLO, FANC, ZMIN, ZMAX 

901 FORMAT (1H 9E14.7) 

C 

C If contour increment has not been set, compute a 'NICE' value 

C for it . 

C 

IF (FANC) 160,170,190 
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]60 CRAT = -FANC 

170 FANC = (HA-GLO ) / ( CRAT-1 ) 

IF (FANC) 220,220,180 

1P0 P = 10.**( IFIX(ALOG10(F/'NC)+500. )-500) 

FANC = AIN'T (FANC/P)*P 

Recompute 'NICE* values of HA and GID. 

GLO = AINT { GDO/ FANC ) * FANC 
HA = A TNT ( HA/ FANC ) * FANC 

Compute contour levels array. 

CO 200 K=1,NLM 
CC = GLO + FLOAT ( K -1 ) * FANC 
KK = K 

IF (CC.GT.HA) GO TO 710 
CL(K) = CC 
200 CONTINUE 
210 NCI, = KK - 1 

Shave away contour values not strictly between ZT1IN and ZMAX. 

230 CONTINUE 

IF (NCL.LE.0) GO TO 240 
IF (CL(1) .GT.ZMIN) GO TO 270 
IF (NCL.LE.l) GO TO 260 
DO 250 1=2, NCL 
250 CI.(I-I) = CL(I) 

260 NCL = NCL - 1 
GO TO 230 
270 CONTINUE 

IF (CL(NCL) .GE.ZMAX) GO TO 260 

240 CONTINUE 
NVAL = 10 

IF (NCL.LT.10) NVAL. = NCL 
DO 741 1=1, NVAL 
J = NCL - I + 1 
CVAL(I) = CL{J ) 

241 CONTINUE 
RETURN 

77.0 ICNST = 1 
NVAL = 10 

IF (NCL.LT.10) NVAL = NCL 
DO 271 1=1 , NVAL 
J = NCL - I •+ 1 
CVAL(I) = CL( J ) 

221 CONTINUE 
RETURN 
C 

END 
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PROGRAM FPLCT2 

DIMENSION X(33) ,Y(33) ,Z(1099) ,M<2178) , 

S (6 ) ,PTA(40) 

DIMENSION MXS (2) ,MXF(2) ,MXJ(?) ,MYS(2) ,MYF(2) ,MYJ( 


901 

900 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


* 

★ 

* 

* 


* 

★ 


COMMON /SRFI3LK/ L1MU { 1024 ) ,LTML( 1024) ,CL(41) ,NCL 

LL .FACT , 1ROT ,NDRZ 

NUPPER .NRSWT ,BIGD ,UMIN 


UMAX ,VMIN , VMAX , RZERO 

NOFFP , NSFVAIi , SPV .BIGEST 

COMMON /PWRZ1/ XXMIN .XXMAX , YYMIN .YYMAX 
ZZMIN , ZZMAX , DELCRT , EYEX 
EYEY , EYEZ 

COMMON / SET/ CVAL (10) .NCHAR , IHEAD(IO) ,XDIST,NON 


,NVAL 


* , X1G , X2G , Y1G, Y2G , ZCMAX, S . NY, AI J , NX 
COMMON /SETN/ NN 
DATA SPVAL, IOFFP/O.O.O/ 

DATA STEREO /O.O/ 

DATA IFR, ISTP, IROTS, IDRX, IDRY, IDRZ, IUPPER, ISKIRT.NCLA/ 

* 1 , 0 , 0 , ], 1 , 1 , 0 , 0 , 0 / 
DATA THETA, HSKIRT, CHI, (TO, CINC/ 

* . 02 , 0 ., 0 ., 0 ., 0 ./ 

DATA S/-30 . , 0. , 3. , 0. , 0. , 0. /MDZ , NX,NY,NCHAR/33, 33, 33, 10/ 

NN = 1 

OPEN (UNIT=4,NAME=' SCRT.DAT' ,TYPE='OLD' , FORM= ’ UNFORMATTED ' , 
* ACCESS^' DIRECT' , ASSOCIATEVARIABLE^NN, RECORDS IZE=1090) 

READ ( 4 ' NN ) I57EAD, X , Y, XXM IN , XXMAX, YYM IN , YYMAX , ZZM IN , ZZMAX, 


* DELCRT, EYEX, EYEY, EYEZ ,NCL, LL, FACT, IROT, NDRZ , NUPPER, NRSWT, 

* BIGD, UMIN, UMAX, VMIN, VMAX, RZERO, NOFFD, NSPVAL, SPV, BIGEST, CL 


* , NX, NY, MDZ , XI G , X2G , YIG , Y2G , ZCMAX, AI J 
READ {4 ’NN){LIMU(l) ,1=1 , 1024) 

READ (4 'NN ) {LTML( I ) , 1=1 , 1024 ) 

READ (4‘NN)(M(I), 1=1, 2178) 

READ (4 'NN) ( Z( I ) , 1=1 , 1089 ) 

READ ( 4 ' NN ) NXP ASS , NYP AS S , MXS , MXF , MXJ ,MYS, MYF , MYJ , RX , RY , I DAY , 


* HSKIRT, NDRZ,XDIST, CVAL, NCHAR, NON, NVAL, S,NNXX, NNYY,MMXX, ISKJRT 
CALL SET3? (X(l),Y(l),Z(l),Ur,VT,DU,4) 

CLOSE (UNIT=4) 

FORMAT (1H 9E14.7) 

FORMAT (1H 13110) 

PRINT 901, X,Y 
PRINT 901, Z 

PRINT 901 , XXMIN, XXMAX, YYMIN, YYMAX, ZZMIN, ZZMAX, DELCRT 
PRINT 901, EYEX, EYEY, EYEZ, FACT, BIGD, BIGEST, RZERO 
PRINT 901, UMIN, LMAX, VMIN, VMAX 

PRINT 900, IHEAD, NCL, LL, IROT, NDRZ , NUPPER, NRSWT, NOFFP, NSPVAL 
PRINT 900.LIMU 

print doo.ijml 

PRINT 901, RX.RY, HSKIRT, XDIST, CVAL, S, CL 

PRINT 900 , NXPASS , NYPASS , MXS , MXF , MXJ , MYS , MYF , MYJ , I DAY, NCHAR 
PRINT 900,f01,NVAL, NNXX,NtTYY,MMXX 
GO TO 102 

do si i=i,mxx 

DO 50 J=l, 8 
K=8+J 


L=16+J 
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N=24+J 

C PRINT 903,1, J,X(I) ,Y{J) ,Z(I,J) ,I,K,X(I) ,Y(K) ,Z(I,K) , 

c * i,l,x(i),y(l),z(i,l),i,n,x(i),y{n),z(i,n) 

903 FORMAT (1H 4 (213, 3 (2X,F5. 2) , 5X) ) 

50 CONTINUE 
PRINT 905 

51 CONTINUE 

905 FORMAT (1H ) 

DO 70 1=1 ft'NXX 
DO 75 J=l,f? 

K=fl4J 

L=1G-KJ 

N=24-HJ 

C PRINT 904, 1,J,X(I),Y(J),M(1,I,J),I,K,X(I),Y(K),M{1,I,K), 

C * I , L,X(l ) ,Y(L) ,M(1, 1, L) ,I,N,X(I) ,Y(N) ,M( 1, 1,N) 

904 FORMAT (1HC 4 (213, 2(2X,F5. 2 ) , 15, 5X) ) 

75 CONTINUE 
PRINT 905 

76 CONTINUE 

DO 101 1=1, NNXX 

DO 100 J=l,8 

K=R+J 

L=16-KJ 

N=24-kJ 

C PRINT 904, I,J,X(l),Y(J),M(2,I,J),I,K,X(I),Y(K),M(2,I,K), 

C * I,L,X(l),Y(L),M(2,I,L),I,N,X(l),Y(N),M(2,I,N) 

100 CONTINUE 
PRINT 905 

101 CONTINUE 

102 CONTINUE 

NON = 1 

OPEN (UNIT=3 , NAME= ' PLOT . DAT ' , TYPE= ’ NEW ' , FORM= 1 UNFORMATTED ' , 

* ACCESS=' DIRECT' ,ASSOCIATEVARIABLE=NON, RECORDS IZE=10) 

120 IF (NXPASS.EQ.2.AND.NYPASS.EQ.2) GO TO 146 
IF (ISKIRT.EQ.O) GO TO 126 

IN = MXS(l) 

IF = MXF(l) 

JN = MYS(l) 

JF = MYF(l) 

IF (NYPASS.NE.l) GO TO 123 

CALL SET32 (X ( 1 ) , Y( JN ) , HSKIRT, UX1 , VX1 , DUMMY, 2 ) 

CALL SET32 (X(NNXX) , Y(JN) , HSKIRT, UX2, VX2, DUMMY, 2) 

DO 121 I=1,M9XX 

CALL SET32 (X ( I ) , Y ( JN ) , KSKI RT, UX , VX , DUMMY, 2 ) 

MU = UX 
MV = VX 

GALL DRAN(MU,MV,M(2*I-1+2*NNXX*(JN-1 ) ) ,M(2*I+2*NNXX*(JN-1) ) ,1 ) 

121 CONTINUE 

CALL DRAW ( IFIX(UN1 ) , IFIX(VX1 ) , IFIX(UX2 ) , IFIX{VX2 ) , 3 ) 

IF ( IDRY.NE.O) GO TO 123 
DO 122 1=2, NNXX 

CALL DRAW{ M( 2*I-3+2*NNXX* ( JN-1 ) ) ,M(2*I-2+2*NNXX*(JN-l ) ) , 

* M( 2*1 -1+2*NNXX*( JN-1 ) ) ,M{2*I+2*NNXX*(JN-1 ) ) ,3) 

122 CONTINUE 

123 IF (NXPASS.NE.l) GO TO 126 
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CALL SET32(X(IN),Y(1),ISKIRT,UY1,VY1,DUMMY,2) 

CALL SE7T32 (X(IN) , Y(NNYY) ,HSKIRT f UY2,VY2, DUMMY, 2) 

DO 124 J=1 , MIYY 

CALL SET32 (X{IN) ,Y(J) , HSKIRT, UY, VY, DUMMY, 2) 

MU = UY 
MV = VY 

CALL DRAW(MU,MV,M(2*IN-1+2*NNXX*(J-1 ) ) ,M(2*IN+2*NNXX*(J-1) ) , 1 ) 

124 CONTINUE 

CALL DRAW( IFIX(UY1 ) , IFIX(VY1 ) , IFIX(UY2 ) , IFIX (VY2 ) , 3 ) 

IF (IDRX.NE.O) GO TO 126 
DO 125 J = 2.NNYY 

CALL DRAJ7(M(2*IN-1+2*NNXX*(J“2)) ,M(2*IN+2*NNXX*(J-2) ) , 

* M{2*IN-1+2*NNXX*(J-1)) ,M(2*IN+2*NNXX*(J-1) ) ,3) 

125 CONTINUE 

126 LI = MXJ(l) 

MI = MXS(l) - LI 
NI = IABS(MI-MXF { 1 ) ) 

LI = MYJ(l) 

MJ = MYS(l) - LJ 
NJ = IABS{MJ-MYF(1)) 

IF (ABS(FX) .LE.ABS(RY) ) GO TO 133 
IF (ISKLRT.NE.O.OR.NYPASS.NE.1) GO TO 128 
I = MXS(l) 

DO 127 J=2,M1YY 

CALL DRAW(M(2*T-l+2*NNXX*(J-2) ) / M(2*I+2*NNXX*{ J-2 ) ) , 

* M(2*I— 1+2*NNXX*(J-1 ) ) ,M(2*I+2*NNXX*(J-1) ) ,2) 

127 CONTINUE 

128 DO 132 II = 1 ,NNXX 
I = MI + II*LI 
IPLI = I + LI 

IF (NYPASS.EO.l) GO TO 129 
K = MYS(l) 

L = MYS (2 ) 

C PRINT 900,NDRZ, II, NI 

IF (IDRX.NE.O) 

* CALL DRAU(M(2*I-1+2*NNXX*(K-1 ) ) ,M(2*I+2*NNXX*(K-1) ) , 

* M(2*I-1+2*NNXX*(L-1)) ,M(2*I+2*NNXX*(L-1 ) ) , 3) 

IF (NDRZ .NE. O.AND. II .NE.NI) 

* CALL CTCELL(Z,MMXX,NNXX,NNYY,M,MINO(l,LfLl),K) 

129 DO 131 JPASS=1 , NYPASS 

LJ = MYJ(JPASS) 

MI = MYS(JPASS) - LJ 
NJ = IABS(MI-MYF{ JPASS) ) 

DO 130 JJ = 1 , NJ 

J = MJ + JJ*LJ 
JPLJ = J + LJ 

IF (IDRX.NE.O. AND. JJ.NE.NJ) 

* CALL DRAW(M{ 2*I-1+2*NNXX*{ J-l ) ) ,M( 2*1+2 *NNXX* ( J-l ) ) , 

* M( 2*I-1+2*NNXX* ( JPLJ-1 ) ) ,M ( 2*I+2*NNXX* ( JPLJ-1 ) ) , 3 ) 

IF (i.NE.MXF(l) .AND. IDRY.NE. 0) 

* CALL DRAW(M(2*IPLI-1+2*NNXX*(J-1)),M(2*IPLI+2*NNXX*(J-1)), 

* M(2*I-1+2*NNXX*(J-1 ) ) ,M(2.*I+2*NNXX*( J-l) ) , 3) 

C PRINT 900,NDRZ,JJ,NJ,II,NNXX,I,J,LI,LJ 

IF(NDRZ.NE. O.AND. JJ.NE.NJ. AND. II 

* .NE.NNXX) 
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* CALL CTCELL(Z,MMXX # NNXX,NNYY,M, 

* MINC(I, I+LI) ,MINO(J, J+LJ) ) 

130 CONTINUE 

131 CONTINUE 

132 CONTINUE 
GO TO 140 

133 IF (ISKIRT.NE. 0. OR. NXPASS. NE. 1 ) GO TO 135 
J = MYS ( 1 ) 

DO 134 I=2,NNXX 

CML DRAW(M(2*I— 3+2*NNXX*(J-l ) ) ,M{2*I-2+2*NNXX*(J-l ) ) , 

* M( 2*1-1 +2*NNXX*( J-l ) ) ,M(2*I+2*NNXX*(J-1 )) ,?) 

134 CONTINUE 

135 DO 13? JJ = 1 ,NNYY 

J = MJ + JJ*LJ 
JPLJ ~ J + LJ 
IF (NXPASn.EQ.l) GO TO 136 
K - MXS(l) 

L = MXS(2) 

IF (IDRY.NE.O) 

* GALL DRAW(M(2*K-1+2*NNXX*(J~1 ) ) ,M(2*K+2*NNXX*(J-1 ) ) , 

* M( 2*L-1+2*NNXX*( J-l ) ) ,M(2*L+2*NNXX* ( J-l ) ) , 2) 

IF (NDRZ .NF. 0. AND. JJ .NE.NJ ) 

* CALL CTCELL ( Z , MMXX , NNXX, NNYY, tt , K, MINO ( J , J+LJ ) ) 

136 DO 130 IPASS=1 , NXPASS 

LI = MXJ(IPASS) 

MI = MXS ( I PASS ) -LI 
NT = IADS (MI -MXF ( I PASS ) ) 

DC 137 11=1, NI 

I = 111 + II*LI 
I PL I = I + LI 

IF( IDRY.NE.O. AND. II. NE.NJ ) 

* CALL DRAW(M(2*I-1+2*NNXX*(J-1) ) ,M(2*I+2*NNXX*(J-1) ) , 

* K( 2*IPLI-1+2*NNXX*( J-l ) ) ,M( 2*IPLI+?*NNXX* ( J— 1 ) ) , 3 ) 

IF (J.NE.MYF(l) .AND.IDRX.NE.O) 

* CALL DRAW(M( 2*1-1 +2*NNXX*(JPLJ-1 ) ) ,M( 2*1+2 *NNXX*(JPLJ-1 ) ) , 

* M(2*I-1+2*NNXX* (J-l ) ) ,M( 2*I+2*NNXX*(J— 1 ) ) , 3) 

C PRINT 900 , NDR7. , II , NI , JJ , NNYY 

IF(NDRZ.NE.O.AND.II.NE.NI.AND.JJ.NE.NNYY) 

* CALL CTCELL NNYY, M, 

* MIN0(I, I+LI) ,MIN0(J, J+LJ)) 

137 CONTINUE 

13P CONTINUE 

139 CONTINUE 

140 IF (ISKIRT.FO.O) GO TO 149 
IF (IDRX.NE.O) GO TO 143 
DO 142 I PASS=1, NXPASS 

IF (NXPASS. EQ- 2) IF=l+( I PASS -1 )* (NNXX-1 ) 

DO 141 J=2,M4YY 

CALL PRAM(M(2*IF-l+2*NNXX*(J-2) ) ,M(2*IF+2*NNXX*(J-2) ) , 

* M(2*IF-1+2*NNXX*(J-1 ) ) ,M(2*IF+2*NNXX*(J-1 ) ) ,1) 


141 


CONTINUE 

142 

CONTINUE 

143 

IF 

(IDRY.NE.O) GO TO 1^9 


DO 

145 JPASS=1 , NYPASS 


IF (NYPASS . EO • 2 ) JF=1+ (JPASS-1 ) * (NNYY-1 ) 
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DO M4 I=2,NNXX 

CALL DRAWM2*I-3+2*NNXX*(JF-l)),M(2*I~2+2*NNXX*(JF-l)), 

* M( 2*1-1 +2*NNXX*(JF-1 ) ) ,M(2*I+2*NNXX*(JF-1 ) ) , 1 ) 

144 CCNTINUE 

145 CONTINUE 
GO TO 149 

146 IF (NUPPER. GT. 0. AND. S (3 ) .LE.S{6 ) ) GOTO 149 
IF (NUPPER. LE. O-AND. S(3) .GT. S(6)) GOTO 149 
NUPPER = 1 

IF (f.(3).LE.P(6)) NUPPER = -1 
DO 14f< I=1,M\D<X 

DO 147 J=1 , NNYY 

IF ( IDRX.NE.C. AND. J.NE. NNYY) 

* CALL DRAW{M( 2*1-1 +2*NNXX*( J-l ) ) ,M(2*I+2*NNXX*(J-1 ) ) , 

* M( 2*1-1 +2*NNXX*J), M( 2*I+2*NNXX*J ), 1 ) 

IF ( IDRY.NE.C. AND. I .NE.NNXX) 

* G ALL DRAW(M(2*I-1+2*NNXX*(J-1 ) ) f M(2*I+2*NNXX*(J-l) ) , 

* M{ 2*I+1+2*NNXX* ( J-l ) ) ,M( 2*I+2+2*NNXX* ( J-l ) ) , 1 ) 

IF( IDR2.NE. O.AND. I.NE.NNXX.AND. J.NE-NNYY) 

* CALL CTCELL ( Z , MMXX , NNXX, NNYY , M , I , J ) 

147 CONTINUE 

148 CONTINUE 

149 IF(STER.EQ.O. ) GO TO 153 
IF ( ISTP) 151,150,152 

150 CONTINUE 

151 CONTINUE 
GO TO 154 

152 IF(IPIC.NE. 2) GO TO 154 

153 CONTINUE 

154 CONTINUE 

WRITE (3 ’NON) END, IND, I END 
CLOSE (UNIT=3) 

RETURN 

END 
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SUBROUTINE SEG(IX1, IY1, 1X2, IY2,KFIAG) 


C 

C 

C 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

5 

c 

999 


6 


66 


C 


C 

998 

77 

7 


TV) is routine v»ou]d normally be the interface between the SRFACE 
package and the plot package. Because of space limitations in 
our machine it is used to write a file containing the data to 
be plotted. Hie PLOTER program reads this file and plots the 
data . 


LOGICAL* 1 START, TITLE, CCNTUR 

COMMON /SET/CVAL(10) , NCHAR ,IHEAD{10) ,XDIST 

,NON , NVAL, X1G, X2G, Y1G, Y2G , ZCMAX, S, NY, AI J, NX 
DATA START/ . TRUE . / , TITLE/ . TRUE . / , UPERIN/102 . / , XIAST/O . / 
DATA CCNTUR/ . TRUE . / 

INTEGER CHEAD(48) 


* 

* 

* 

* 


'O ' 

, 'TH' 

,'E * 

, ‘FI ‘ , 

' RS ‘ 

,‘T ‘ 

1 1 
9 

1 1 
r 9 

'ON' 

, 'TO' 

, ’UR' 

,' V‘, 

'AL' 

,‘UE’ 

,'S • 

, 'OF', 

'TA' 

, 'NT' 

z' 

,’ (*, 

‘HI* 

, ’GH’ 

, 'ES' 

» rn • 

/ -L 9 

■ow 1 

, 'ES 1 

,'T)' 

i i 

i 9 

' WA ' 

, ' TT ' 

,’S/' 

, ' CM* , 


' ** 1 , ' 2 7 

DIMENSION CVALQ(P) ,S(6) 


On the first time through this routine the program writes all 
the headers and labes for the graph. Daring subsequent calls the 
the hidden line data is written for ploting. 

IF (.NOT. START) GO TO 10 

TYPE 999, NY, X1G,X2G, Y1G, Y2G, ZCMAX 

FORMAT (HI 15, 5F12.4) 

XLAST = 0. 

IXTO = 0 
IYTO = 0 

IF ( .NOT.CCNTUR) GO TO 8 

HT = 8.0 

DO 6 1=1,4 

HT = HT - .5 

II = 12*(I-1)+1 

IK = II + 12 

WRITE (3 'NCN)HT, (CHEAD(J) ,J=II, IK) ,NVAL 
CONTINUE 

WRITE (3 'NON) X1G, X2G, Y1G, Y2G, ZCMAX 
WRITE (3 'NON) S (1 ) ,S(2) ,S(3) ,NY,NX 
HT = 5.5 
DO 7 1=1, NVAL 
HT = HT - .5 
TYPE 999, I,CVAL(I ) ,AIJ 
CVAL(I) = CVAL(I)*AIJ 
ENCODE (12, 77, CVALQ) CVAL( I) 

TYPE 998, CVALQ 
FORMAT (1H PA2) 

FORMAT (FI 2. 1 ) 

WRITE ( 3' NON )HT, CVALQ 
CONTINUE 
CONTUR = .FALSE. 
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P CONTINUE 

IF ( .NOT. TITLE) GO TO 10 
WRITE ( 3 'NON )XDIST> IHEAD 
TITLE = .FALSE. 

10 IF ( IXI .EQ. IXTO.AND. IYI .EQ. IYTO) GO TO 20 
AX = FLOAT( IXI ) /UPERIN 
AY = FLOAT (IYI) /UPERIN 
IP3 = +3 

WRITE (3 ’NON) AX, AY, IP3 
IF (KFLAG.GT.l) PRINT 901, IPC, AX, AY 
901 FORMAT (1H I10,2E14.7) 

20 BX = FLOAT (1X2) /UPERIN 

BY - FLOAT ( IY2 ) /UPERIN 
IP2 = +2 

WRITE (3’N0N)EX,BY, IP2 

IF (KFLAG.GT.l) PRINT 901, IP2,BX,BY 

IXTO = 1X2 

IYTO = IY2 

START = .FALSE. 

XLAST = AMAX1 (XLAST, AMAXO ( 1X1 , 1X2 ) ) 

RETURN 

END 
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c 
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SUBROUTINE SET32 (X,Y,Z, XT, YT,ZT, KFLAG) 

Ibis routine implements the 3 -space to 2-space transformation 
by Kuber, Szabo, and Giulieri, The Perspective Representation 
of Functions of TVo Variables, J. ACM 15, 2 193-204,1968 

X, Y,Z Aie the 3-space coordinates of the intersection of the 

line of sight and the image plane. This point can be 
thought of as the point looked at. 

XT, YT, ZT Are the 3-space coordinates of the eye position. 


C KFLAG = 


2 argunents 


C 

C X, Y, 
C 

C XT, 

C ZT 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Z Are the 3-spoce coordinated of a point to be 

transformed . 

YT The results of the 3-space to 2-space transformation. 
Not used. 

If LL{in OCMMQN)=0, XT and YT are in the same scale and X,Y 

The variable KFLAG has tvvo possible variables 

1- cctnpute intersection of line of sight 

2- transform from 3-space to 2-space 

NOTE! ! 1! !!!!!!! ! 

The KFLAG=3,4 are special debugging flags and are not 
part of the plot package. 


* 

* 

* 

* 

* 

★ 


COMMON 

/PWRZl/ XXMIN 

,XXMAX 

, YYMIN , YYMAX , 


ZZMIN 

, ZZMAX 

, DELCFT , EYEX 


EYEY 

, EYEZ 


CCMMCN 

/SRFBLK/ MMU(1024) 

,LIML(1024) ,CL(41 ) ,NCL 


LL 

, FACT 

, IROT , NDRZ 


NUPPER 

, NRSWT , BIGD ,UMIN 


UMAX 

, VMIN 

, VMAX , RZERO 


NOFFP 

, NSFVAL , SPV , BIGEST 

COMMON 

/SETN/ NN 



DIMENSION NLU(7 ) 

,NRU{7) 

, NB V ( 7 ) ,NTV{7 ) 


c 


C picture comer coordinates for 11^=1 

C 

DATA NLU(1 ) ,NRU(1 ) ,NBV(1) ,NTV(1 )/l0, 1014, 10, 1014/ 


Z. 


picture corner coordinates for LL=2 

DATA NUU{2) ,NRU(2) ,NBV(2) ,NTV(2)/10, 924,50, 964/ 


picture comer coordinates for LL=3 

DATA NUJ(3 ) ,NRU{3 ) ,NBV{3 ) ,NTV(3 )/l00, 1014, 50, 964/ 


picture comer coordinates for LL=4 

DATA NLU(4 ) ,NRU(4) ,NBV(4) ,NTV(4 )/l0, 1014, 10, 1014/ 
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c 

c 

c 

c 

c 

c 

c 

c 

c 


1 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

900 

901 


c 

c 

c 


c 

c 

c 


picture comer coordinates for LL=5 

DATA NLU ( 5 ) , NRU ( 5 ) ,NBV(5 ) ,NTV(5 )/10, 1014, 10, 1014/ 

picture comer coordinates for LL=6 

DATA NIJU{6) ,NRU(6) ,NBV(6) ,OTV{6)/10, 512,256,759/ 

picture comer coordinates for LL=7 

DATA NLU{7) ,NRU(7) ,M3V(7) ,17TV(7 )/512, 1014, 256, 750/ 

GO TO (1,2, 3,4) KFLAG 
JUMP3 =104 

IF (N0FFP.EQ.1) JUMP3 = 103 

AX = X 

AY = Y 

AZ = Z 

EX = XT 

EY = YT 

EZ = ZT 

As much computation as possible is done during execution of 
SET32 since the transformation is called many times for each 
call to SET32. 


DX = AX - EX 
DY = AY - EY 
DZ = AZ - EZ 

A more careful computation of direction Cosines. 


D = 0. 

T = AMAX1 (ABS(DX) ,ABS(DY) ,ABS(DZ) ) 

PRINT 901 , AX, AY, AZ , EX, EY, EZ , DX, DZ, DZ , T 
PRINT 900,NOFFP 
FORMAT (III 120010) 

FORMAT (1H 9E14.7) 

IF (T.EQ.0.0) GO TO 30 
R1 = DX/T 

R2 = DY/T 
R3 = DZ/T 

D = SQRT(R1*R1 + R2*R2 + R3*R3) 

If D isn't ZERO 

C0SAL = Rl/D 
COSBE = R2/D 
COSGA = R3/D 
D = D*T 
GO TO 40 

If D is ZERO, ray has no direction: assign direction down 
X-axis . 
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30 CONTINUE 

COSAL = 1. 

COSBE = 0. 

CCGGA = 0. 

40 CafTINUE 

AL = ACOS( COSAL) 

BE = ACOS (COSBE) 

GA = ACOS ( COSGA ) 

SINGA = SIN(GA) 

C PRINT 901, R1,R2,R3,D, COSAL, COSBE, COSGA,AL, BE, GA, SINGA 

C PRINT 900, LL 

JUMP2 = 110 
IF (LL.EQ.O) GO TO 101 
JUMP2= 103 

DELCRT = NRU(LL)-NLU(LL) 

UO = UMIN 
VO = VMIN 
U1 = NLU(LL) 

VI = NBV(LL) 

U2 = NRU(LL) - NLU(LL) 

V2 = NTV(LL) - NEV(LL) 

IF(UMAX-UMIN) 52,51,52 

51 U3 = 0. 

GO TO 53 

52 U3 = U2/(UMAX-UMIN) 

53 CONTINUE 

IF (VMAX-VMIN) 55,54,55 

54 V3 = 0. 

GO TO 56 

55 V3 - V2/ (VMAX-VMIN) 

56 CONTINUE 

U4 = NRU(LL) 

V4 = NTV(LL) 

C PRINT 901,UO,VO,U1,V1,U2,V2,U3,V3,U4,V4,LMAX,VMAX 

C PRINT 900,NRSWT 

IF (NRSWT.BQ.O) GO TO 101 
UO = -BIGD 
VO = -BIGD 
U3 = U2/(2.*BIGD) 

V3 = V2/(2.*BIGD) 

C 

C Tie 3-space point looked at is transformed inot (0,0) of the 

C 2-space. The 3-space Z axis is transformed into the 2-space Y 

C axis. If the line of sight is close to parallel to the 3-space 

C Z axis, the 3-space Y axis is chosen (instead of the 3-space 

C Z axis) to be transformed into the 2-space Y axix. 

C 

101 IF (SINGA.LT. 0.0001) GO TO 102 
R = 1. /SINGA 
JUMP = 105 

C PRINT 900, JUMP, JUMP2,JUMP3 

RETURN 

102 SINBE = SIN (BE) 

R = 1. /SINBE 
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JUMP = 106 

C PRINT 900, JUMP, JUMP3,JUMP2 

RETURN 
C 

C Transformation entry point 

C 

2 XX = X 
YY = Y 
ZZ = Z 

C PRINT 901,XX, YY, ZZ 

C PRINT 900, JUMP3 

IF (JUMP3.EQ.104) GO TO 104 

103 IF (Z2.EQ.SPV) GO TO 109 

104 DENCH = (XX -EX ) *OCSAL + ( YY-EY) *COSBE + ( ZZ-EZ )*C0SGA 
IF (DENCM.I3E.0.0) GO TO 1111 

Q - 1. 

GO TO 50 

1111 Q = D/DEN0M 
50 CONTINUE 
C PRINT 901,DENCM,Q 

C PRINT 900, JUMP 

IF ( JUMP. EQ. 106) GO TO 106 

105 XX = ( (EX4Q*(XX-EX)-AX)*00SBE-(EY4Q* (YY-EY) -AY)* COSAL )*R 
YY = ( EZ+Q* ( ZZ-EZ ) -AZ ) *R 

GO TO 107 

106 XX = ( (EZ4Q*( ZZ-EZ ) -AZ )* COSAL- ( EX4Q*(XX-EX) -AX) *COSGA)*R 
YY = ( EY4Q* (YY-EY ) -AY ) *R 

107 IF (JUMP2.EQ.110) GO TO 110 

108 XX = AM INI (U4, AMAX1 (Ul , Ul+U3*( FACT*XX-UO) ) ) 

YY = AMIN 1 ( V4 , AMAX1 (VI , V1+V3 * ( FACT*YY-VO ) ) ) 

GO TO 110 

109 XX = NSPVAL 
YY = NSPVAL 
GO TO 110 

no xt = xx 

YT = YY 

C PRINT 901 ,XT, YT 

C PRINT 900, JUMP2 

RETURN 

3 CCNTINUE 

WRITE (3 'NN ) JUMP, JUMP2, JUMP3, EX, EY, EZ, COSAL, COSBE, COSGA, D, 

* AX,AY,AZ,R,UO,U1,U2,U3,U4,VO,V1,V2,V3,V4 

C PRINT 900, JUMP, JUMP2,JUMP3 

C PRINT 901 , IX, EY, EZ , COSAL, COSEE, COSGA , D, AX, AY, AZ , R 

C PRINT 901, UO, U1,U2,U3,U4, VO,Vl, V2,V3,V4 

RETURN 

4 CONTINUE 

READ (4 'NN) JUMP, JUMP2,JUMP3, EX, EY,EZ, COSAL, COSDE, COSGA, D, 

* AX, AY, AZ, R,U0, Ul ,U2, U3, U4, VO, VI , V2, V3, V4 

C PRINT 900, JUMP, JUMP2.JUMP3 

C PRINT 901, EX, EY.EZ, COSAL, COSBE, COSGA, D,AX,AY,AZ,R 

C PRINT 901,UO,U1,U2,U3,U4,VO,V1,V2,V3,V4 

RETURN 
END 

FLNCTICN ACOS(X) 
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ACOS = 1 . / SQRT ( 1 . +ATAN (X ) * ATAN (X ) ) 
END 

FUNCTION R(HO,HU) 


TViis routine interpolates in the CV array. 


COMMON /CVAL/CV 
IF (IKMIU.EQ.O) GO TO 10 
R ~ ( HO-CV ) / ( HO— HU ) 

C PRINT 900, CV, HO, HU, R 

900 FORMAT (1H 4(2X,E14.7)) 
RETURN 

10 CONTINUE 

IF (HO-CV. LT. 0.0) R = 0.0 
IF (HO-CV. GE. 0.0) R = 1.0 
RETURN 

era 
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c 

c 

c 

c 

c 

c 

c 

c 

c 


1 

2 

3 

4 


C 

901 


101 

c 

c 

c 

c 

c 


c 


c 


c 

c 

904 

905 
10 


SUBROUTINE CTCELL (Z,MDZ,NX,NY,M, 10, JO) 

CTCELL ccnputes lines for constant Z (coutour lines) in one 
cell of the array Z for the SRFACE package. 

Z, 1]X, MY are the same as in SRFACE. 

By the time ctcell is first called, M contains 
the two-space plotter location of each Z point . 
U(Z(I,J))=M(1,I,J), V(Z(l,J))=M(2,I,J) 

10, JO Hie cell Z(I1, Jl) to Z(I1+1,J1+1) is the one to 
be contoured . 

DIMENSION Z( 1089) ,M(217P ) 

COMMON /SRFBLK/ LIMU(1024) ,LIML{1024) ,CL{41) ,NCL 
LL ,FACT , IROT , NDRZ 

NUPPER , NRSOT , BIGD ,UMIN 

UMAX , VMIN ,VMAX , RZERO , 

NOFFP , NSPVAL , SPV , BIGEST 

COMMON /CVAL/ CV 
II = IO 
I1P1 =11+1 
Jl = JO 
J1P1 = Jl + 1 
HI = Z(ll-fNX*(Jl-l) ) 

H2 = Z(I14NX*(J1P1-1)) 

113 = Z(I1P14NX*(J1P1-1)) 

H4 = Z( I 1P14NX*(J1— 1 ) ) 

PRINT 901, HI, H2, JI3, H4 
FORMAT {1H 9E14.7) 

IF (N0FFP.NE. 1 ) GO TO 101 

IF ( HI . EQ. SPV . OR. H7 . EQ. SPV.OR .H3 . EQ. SPV. OR. H4 . EQ. SPV ) RETURN 
IF (TMIN1 (HI, EI2, H3, H4) .GT.CLfNCL) ) RETURN 

For each contour level, decide vhich of the 16 basic 
siturations exists, then interpolate in two-space to find 
the end points of the contour line segment within this cell. 

DO 111 K=1 , NCL 
CV = CL(K) 

PRINT 901, CV 

K1 = (lFIX(SIGN(l.,Hl-CV))+l)/2 
K2 = { IFIX(SIGN{1 . , H2-CV) )+l )/2 
K3 = {IFIX(SIGN{1. ,H3-CV))+l)/2 
K4 = { IFIX(SIGN (1 . , H4-CV) )+l )/2 
JUMP = 1+K1+K2*2+K3*4+K4*8 
PRINT 900, K1,K2,K3,K4, JUMP 
I FUG = 4 

IF ( JUMP. EQ.1.0R. JUMP. EQ. 16) GO TO 10 
PRINT 904, JUMP, 11, I1P1, Jl, J1P1 
PRINT 905, K1 , HI , K2, H2, K3, H3,K4, H4.CV 
FORMAT (1H 5110) 

FORMAT (1H 4{I3,2X,F5.2) ,2X,F5. 2) 

CONTINUE 

A1 = FIOAT(M(2*Il-l+2*NX*(Jl-l))) 

A2 = FLOAT(M( 2*1 1+2*NX* ( Jl— 1 ) ) ) 

A3 = FLQAT(M( 2*I1-1+2*NX*( J1P1-1 ) ) ) 

A4 = FLDAT(M(2*I1+2*NX*(J1P1-1))) 
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A5 = FIDAT(M(2*I1P1-1+2*NX*(JL-1))) 

A6 = FIDAT(M( 2*IlPl+2*NX*( Jl-1 ) ) ) 

A7 = FIDAT(M(2*I1P1-1+2*NX*(J1P1-1))) 

El = FLOAT (M( 2*1 1-1+2*NX*(J1P1-1 ) )-M( 2 *11-1+2 *NX*( Jl-1 ) ) ) 

B2 = FIOAT(M(2*Il+2*NX*(JlPl— 1 ) )-M(2*Il+2*NX*(Jl-l ) ) ) 

A JR = FLOAT(M(2*IlPl+2*NX*(JlPl-l ) ) ) 

B31 = FIOAT(M( 2*I1P1-1+2*NX*( J1P1-1 ) )-M (2*1 1-1+2 *NX*(J1P1-1 ) ) ) 
B41 = FLOAT(M(2*IlPl+2*NX*( J1P1-1 ) )-M( 2*I1+2*NX*(J1P1-1 ) ) ) 

Bll = FLOAT(M(2*IlPl-l+2*NX*(Jl-l ) )-M(2*Il-l+2*NX*( Jl-1 ) ) ) 

B21 = FL0AT(M(2*I1P1+2*NX*(J1-1 ) )-M(2*Il+2*NX*(Jl-l) ) ) 

B3 = FI0AT(M(2*I1P1-1+2*NX*(J1-1))-M(2*I1P1-1+2*NX*(J1P1-1))) 
B4 = FLOAT (M( 2*1 1P1+2*NX* (Jl-1 ) )-M( 2*I1P1+2*NX*( J1P1-1 ) ) ) 

IF (JUMP.EQ. 1 ) GO TO 112 

IF (JUMP.EQ. 2) GO TO 103 

IF (JUMP.EQ. 3) GO TO 105 

IF (JUMP.EQ. 4) GO TO 106 

IF (JUMP.EQ. 5) GO TO 107 

IF (JUMP.EQ. 6) GO TO 102 

IF (JUMP.EQ. 7) GO TO 108 

IF (JUMP. EQ. 8) GO TO 109 

IF (JUMP.EQ. 9) GO TO 109 

IF (JUMP.EQ. 10) GO TO 108 
IF (JUMP.EQ.il) GO TO 104 

IF (JUMP.EQ. 12) GO TO 107 

IF (JUMP.EQ. 13) GO TO 106 

IF (JUMP.EQ. 14) GO TO 105 

IF (JUMP.EQ. 15) GO TO 103 

IF (JUMP.EQ. 16) GO TO 111 

102 IDUB = 1 

103 RA = R(H1,H2) 

MUA = A1 + RA*B1 
MVA = A2 + RA*B2 
RB = R(H1,H4) 

MUB = A1 + RB*B11 
MVB = A2 + RB*B21 
GO TO 110 

104 IDUB = -1 

105 RA = R(H2, HI) 

MUA = A3 - RA*B1 
MVA = A4 - RA*B2 
RB = R(H2, H3) 

MUB = A3 + RB*B31 
MVB = A4 + RE*B41 
GO TO 110 

106 RA = R(H2,I13) 

MUA = A3 + RA*B31 
MVA = A4 + RA*B41 
RB = R(H1, H4) 

MUB = A1 + RB*B11 
MVB = A2 + RB*B21 
GO TO 110 

107 RA = R(H3,H2) 

MUA = A7 - RA*B31 
MVA = A8 - RA*B41 
RB = R(H3,H4) 
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MUB = A7 + RR*B3 
MVB = AP + RB*E4 
IDUB = 0 
GO TO 110 

1 OB RA = R(H2,H1) 

MUA = A3 - RA*B1 
MVA = A4 - RA*B? 

RB = R(II3,H4) 

MUB = A7 + RB*B3 
MVB - A8 + RE*B4 

GO TC no 

109 RA = R(H4, HI ) 

MUA = A5 - PA* Dll 
MVA = A6 - RA*B21 
RB = R(H4,H3) 

MUB = A5 - RB*B3 
MVB = A6 - RE*B4 
IDUB = 0 

110 CONTINUE 
PRINT 901, RA, RB 

PRINT 900, MUA, MVA., MUB, MVB, IDUB, IFLG 
CALL DRAW (MUA, MVA, MUB, MVB , 1 ) 

IFLG = 4 

900 FORMAT (1H 13110) 

IF (IDUB) 109,111,107 

111 CONTINUE 

112 RETURN 
END 


A-32 



c 

c 

c 

c 
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c 
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* 

* 

* 
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55 

c 

C 

c 

1 

c 

900 

C 

C 

C 

2 


C 

c 

c 

3 

c 

c 

c 

101 


SUBROUTINE DRAW ( MX1 , MY1 , MX? , MY2 , KFLAG ) 

This routine draws the visible part of the line connecting 
(MX1,MY1) and (MX?,MY2). The variable KFIAG is used to 
specify which mode the subroutine uses: 

1-Draw visible part of line 
? -MARKs the visibility arrays 
3-Roth marks and draws 

L0GICAL*1 VIS! , VIS? 

COMMON / SRFBIiK/ LIMU(1024) ,LIML(1024) ,CL(41) ,NCL , 

LL .FACT .TROT ,NDRZ 

NUPPER , NRSWT , BIGD ,IMIN 

UMAX , VT1IN , VMAX .PZERO , 

NOFFP , NSPVAL , SPV , DIGEST 

DATA STEEP/5./ 

ITFG = 1 

IF { KFLAG . LT . 4 } GO TO 55 
ITFG = KFLAG 
KFLAG = KFIAG - 3 
GO TO (1,2,3) KFIAG 

DRAW 

IDRAW = 1 
IMARK = 0 

PRINT 900, MX1 , MY1 , MX? , MY2 , KFIAG, IDRAW, IMARK 
FORMAT (1HC 12110) 

GO TO 101 

MARK 

IDRAW = 0 
IMARK - 1 
GO TO 101 

DRAW and MARK 

IDRAW = 1 
IMARK = 1 


MARK line left to right. 

MMX1 = MX1 
MMY1 = MY1 
MMX2 = MX2 
MMY? = MY? 

IF (MMX1.EQ. NSPVAL. OR. MMX2.EQ. NSPVAL) RETURN 
L0GICAL*1 BAD 
BAD = MX1.GT.1024 
BAD = MY1 . GT . 1024 . OR . BAD 
BAD = MX2.GT.1024.0R.BAD 
BAD = MY2.CT.1024.OR. BAD 
BAD = MX1.LT.1 .OR. BAD 
BAD = MY1.LT.1 .OR. BAD 
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BAD KJC2.LT. 1 .OR. BAD 

BAD = MY2.LT. 1 .OR. BAD 

IF (BAD) RETURN 

IF (M1X1.GT.MMX2) GO TO 102 

NX1 - MMX1 

NY1 = MMY1 

NX7 = MKX2 

NY2 = MMY2 

GO TO 103 

102 NX1 = MMX2 
NY1 = MMY2 
NX2 = MMX1 
NY2 = MMY1 

103 IF (NUPPER.LT. 0) GO TO 110 
C 

C Check upper visibility. 

C 

VIS1 = NYI.GE. (LIMUfNXl )— 1 ) 

VIS2 = NY2.GK. (LIMU(NX? )-l ) 

C PRINT 900,NX1,NY1,NX2,NY2,VIS1,VIS2,LIMU(NX1 ) ,LLMU(NX2) 

C 

C VIS1 and VIS2 TRUE means visible. 

C 

IF (VIS1.A1JD.VIS2) GO TO 113 
C 

C VLSI or VIS2 fa] se means invisible. 

C 

IF ( .N0T.(VIS1.CR.VIS2)) GO TO 119 
C 

C Find change point . 

C 

IF (NX1.EQ.NX2) GO TO 112 

DY = FLOAT(NY2-NYl ) /FLOAT (NX2-NX1 ) 

NX1P1 = NX1 + 1 
FNY1 = NY1 
IF (VIS1) GO TO 107 
DO 104 K=NX1P1,NX2 
MX = K 

MY = FNYl + FLOAT (K-NXl)*DY 
IF (MY.GT.LIMU(K)) GO TO 105 

104 CONTINUE 

105 IF (ABS(DY) .GE. STEEP) GO TO 110 

106 NX1 = MX 
NY1 = MY 
GO TO 113 

107 DO 108 K=NX1P1,NX2 
MX = K 

MY = FNYl + FLOAT ( K-NX1 ) *DY 
IF (MY.LT.LIMU(K) ) GO TO 109 

108 CONTINUE 

109 IF (ABS(DY) .GE. STEEP) GO TO 111 

NX2 = MX 
NY2 = MY 
GO TO 113 

110 IF (LIMU(MX) .EQ.O) GO TO 106 
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non non 


NXl = MX 
NYl =LIMU(NX1) 

GO TO 113 

111 NX2 = MX 

NY2 = LIMU(NX2) 

GO TO 113 

112 IF (VIS1) NY2=MIN0{L1MU{NX1 ) ,LIMU(NX2 ) ) 
IF (VIS1) NY1=MIN0(LIMU(NX1 ) ,LIMU(NX2) ) 

113 IF (IDRAW.EQ.O) GO TO 116 


Draw visible part of line. 


IF (IROT) 114,115,114 

CALL SEG(NY1, 1023-NX1,NY2,1023-NX2, ITFG) 
GO TO 116 

CALL SEG (NXl, NYl, NX2,NY2, ITFG) 

IF (IMARK.EQ.O) GO TO 119 
IF (NX1.EQ.NX2) GO TO IIP 
DY = FLOAT (NY2-NY1) /FLOAT (NX2 -NXl) 

FNY1 = NYl 


DO 117 K=NX1 ,NX2 

LIMU(K) = FNY1 + FLOAT (K-NX1)*DY 

117 CONTINUE 
GO TO 119 

118 LIMU(NXl) = MAXO (NYl , NY2 ) 

IF (NUPPER) 120,120,138 


Same idea as above, but for lower side. 

120 IF (mXl.GT.MMX2) GO TO 121 
NXl = MMX1 

NYl = MMY1 
NX2 = MIX2 
NY2 = mY2 
GO TO 122 

121 NXl = MMX2 
NYl = MMY2 
NX2 = MMXl 
NY2 = MMY1 

122 VTS1 = NY1.LE. (LIML(NXl)tl ) 

VIS 2 = NY2.LE. (LIML(NX2)+1) 

IF (VIS1.AND. VIS2) GO TO 132 

IF ( .N0T.(VIS1.CR.VIS2)) GO TO 138 
IF (NX1.EQ.NX2) GO TO 131 
DY = FLOAT (NY2-NY1 ) /FLOAT (NX2-NX1 ) 
NX1P1 = NXl + 1 
FNY1 = NYl 
IF (VIS1) GO TO 126 
DO 123 K=NX1P1,NX2 
MX = K 

MY = FNY1 + FIOAT(K-NXl )*DY 
IF (MY.LT.LIML(K)) GO TO 124 

123 CONTINUE 

124 IF (ABS(DY) .GE. STEEP) GO TO 129 

125 NXl = MX 
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NY1 = MY 
GO TO 13? 

1 26 DO 127 K=NX1P1,NX2 

MX = K 

MY - FNY1 + FLOAT { K-NXl ) * DY 
IF (MY.Gr.LIML(K) ) GO TO 12R 

1 27 CONTINUE 

128 IF (ABS(DY) .GE. STEEP) GO TO 130 
NX2 = MX 

NY? = MY 
GO TO 132 

12? IF (LIML(MX) .EQ.1024) GOTO 125 
NX1 = MX 
NY1 = LIML(NX1 ) 

GO TO 132 

130 NX2 = MX 

NY2 = LIML(NX2 ) 

GO TO 132 

131 IF (VIS1 ) NY? = MAX0{LIML(NX1) ,LIML(NX2)) 
IF (VIS?) NY1 = MAX0(LIML(NX1),LIML(NX2)) 

132 IF (1DRAU.EQ.0) GO TO 135 
IF (IROT) 133,134,133 

133 CALL SEG(NY1, 1023-NX1 ,NY2, 1023-NX?, ITFC) 
GO TO 135 

134 CALL SEG(NX1,NY1,NX2,NY2, ITPG) 

135 IF (IMARK.EQ.O) GO TO 138 
IF (NX1.EQ.NX2) GO TO 137 

DY = FLOAT(NY2-NY1)/FIOAT(NX2-NX1) 

FNY1 = NY1 
DO 136 K=NX1 ,NX2 

LIML(K) = FNY1 + FLOAT (K-NXl )*DY 

136 CONTINUE 
RETURN 

137 LIML(NXl) = MIN0(NY1,NY2) 

138 RETURN 
END 
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PROGRAM PLOTER 
LOGICAL* 1 FIRST, YY, IY 

INTEGER XIGA(7 ) ,X2GA( 7 ) ,Y1GA(7 ) , Y2GA{ 7 ) ,XAA( 17 ) ,XBA( 17) , 

* XCA(17) ,XDA(17) 

INTEGER CHEAD(48) ,IHEAD(10) ,ALP, AG'I,A\RP 
DIMENSION CVALQ(8),S(3),JIIEAP(20) 

DATA YY/'Y'/ 

DATA ALP/' C/rARP/') '/,Aai/' , 7, FIRST/ .TRUE./ 

NC = 0 
NN = 0 
DLINE = 4. 

FMA.G = .8 

TYPE THE CURRENT PLOT MAGNIFICATION IS .8' 

TYPE DC YOU UANT TO CHANGE PLOT MAGNIFICATION? (Y/N) ' 
ACCEPT 8050, IY 
9050 FORMAT ( A1 ) 

IF (IY.GT.98) IY = IY - 32 
IF (1Y.NE.YY) GO TO 5 

TYPE *, ' ENTER MAGNIFICATION FACTOR -DECIMAL’ 

ACCEPT 906,FMAG 
906 FORMAT (F12.7) 

5 CONTINUE 

CALL PLOTS (0,0,0) 

CALL PLOT (0.0,0.75, -3) 

HT = 8.0 
NO - 1 

OPEN ( UNIT=4 , NAME= ' PLOT . DAT ' , TYPE= ' OLD ' , FORM- 1 UNFORMATTED ' , 

* ACCESS= 1 DIRECT ' , ASSOC IATEVARIABLE=NO , RECORDS IZ E=1 0 ) 

DO 6 K=l,4 

READ ( 4 ' NO ) HT, (CHEA-D(l) , 1=1, 13) ,NVAL 
C PRINT 900, NVAL 

C PRINT 901, HT 

CALL SYMBOL (0. 0, HT, . I , CHEAD(1 ), 0. 0, 24 ) 

6 CONTINUE 

READ (4 ’NO) X1G.X2G, Y1G, Y2G, ZCMAX 
READ (4'NO) S(1),S(2),S(3) ,NY, NX 
NY = NY - I 
NX ~ NX - 1 
66 DO 7 1=1, NVAL 

READ (4'NO) HT.CVALQ 
C PRINT 900, NVAL 

C PRINT 901, HT 

CALL SYMBOL (0.0, ITT, . 1,CVAK>, 0.0, 12 ) 

7 CONTINUE 

READ (4’NO) XDIST, IHEAD 

TYPE *, ' ENTER 40 CHARACTER PLOT TITLE' 

ACCEPT 9055, JHEAD 
9055 FORMAT (20A2) 

C PRINT 901, XDIST 

CALL SYMBOL (XDIST+1 . , 6 . 5 , . 1 , JHEAD, 0.0,40) 

10 CONTINUE 

AXO = AX 

AYO - AY 

READ (4'NO) AX, AY, IPC 
AX = AX*FMAG + 1 . 
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AY = AY*FMAG - .75 
IF (AXO. NE. 0.0) GO TO 11 
XD = AX 
YD = AY 

11 YCMAX = YMAX 

YMAX = AMAX1 (YCMAX, AY) 

IF (YCMAX. EQ.XMAX) GO TO 12 
XO = AX 

12 IF (NN.EQ.O. AND. FIRST) GO TO 18 

IF (ABS(AXO-AX) .GT.DLINE) NN = NN + 1 
IF (NN.NE.O) GO TO 13 
XB = AX 
YB = AY 

13 IF { ABS ( S ( 1 ) ) . LT . ABS ( S ( 2 ) ) ) GO TO 131 
IF (NN.NE.NX) GO TO 18 

GO TO 132 

131 IF (NNf.NF.NY) GO TO 18 

132 IF (NC.NE.O) GO TO 14 
XC = AX 

YC = AY 
NC = NC + 1 

14 CONTINUE 
XA = AXO 
YA = AXO 

1R FIRST = .FALSE. 

900 FORMAT (1H 13010) 

C PRINT 903 , AX, AY, I PC , NO, NC, AXO, AYO 

903 FORMAT (1H 2F12.4, 317, 2F12.4) 

901 FORMAT (1H 9F12.4) 

IF (IPC.EQ.O) GO TO 20 
IF (AX. EQ- -9999. ) GO TO 20 
CALL PLOT (AX. AY, IPC) 

GO TO 10 

902 FORMAT (F6.3) 

20 CONTINUE 

IF (XA.LT.XC) GO TO 21 

XS = XC 

XC = XA 

XA = XS 

YS = YC 

XC = YA 

YA = YS 

21 CONTINUE 

IF (XD.GT.XB) GO TO 22 

XS = XB 

XB = XD 

XD = XS 

YS = YB 

YB = YD 

YD = YS 

22 CONTINUE 

C TYPE 901 , XA, YA,XB, YB, XC, YC, XD, YD, XlG, X2G, YlG, Y2G 

ENCODE (6,902, XiaA )X1G 
ENCODE (6,902, X2GA )X2G 
ENCODE (6,902, Y1GA ) YlG 
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mono 


ENCODE (6,902, Y2GA ) Y2G 
TYPE 905,X1GA 
TYPE 905.X2GA 
TYPE 905, Y1GA 
TYPE 905, Y2GA 

905 FORMAT (1H 28A2) 

IF (ABS(G (1 ) ) .LT.ABS(S (2 ) ) ) GO TO 200 
IF (S(l).LT.O.O) GO TO 100 
XAA(l) = ALP 
DO 25 1=1,3 

25 XAA(I+1) = X1GA(I) 

XAA(5) = ACM 
DO 30 1=1,3 

30 XAA ( 1+5 ) = Y1GA(I) 

XAA{9) = ARP 
C TYPE 905, XAA 

XBA(l) = ALP 
DO 35 1=1,3 

35 XBA( 1+1 ) = X2GA( I ) 

XBA(5) = ACM 
DO 40 1=1,3 

40 XBA( 1+5 ) = Y1GA(I) 

XBA(9 ) = ARP 
C TYPE 905, XBA 

XCA(l) = ALP 
DO 45 1=1,3 

45 XCA( 1+1 ) = X1GA(I) 

XCA( 5 ) = ACM 
DO 50 1=1,3 

50 XCA(I+5 ) = Y2GA(I) 

XCA(9) = ARP 
C TYPE 905, XCA 

XDA(l) = ALP 
DO 55 1=1,3 

55 XDA( 1+1 ) = X2GA( I ) 

XDA( 5 ) = ACM 
DO 60 1=1,3 

60 XDA(l+5) = Y2GA(I) 

XDA(9) = ARP 
GO TO 500 

100 XAA(l) = ALP 

DO 125 1=1,3 

125 XAA( 1+1 ) = X2GA( I ) 

XAA(5) = ACM 
DO 130 1=1,3 

130 XAA( 1+5 ) = Y2GA(I) 

XAA( 9 ) = ARP 
XBA(l) = ALP 
DO 135 1=1,3 

135 XRA( 1+1 ) = X1GA(I) 

XRA(5) = ACM 
DO 140 1=1,3 

140 XBA{ 1+5 ) = Y2GA(I) 

XBA(9) = ARP 
XCA(l) = ALP 
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DO 145 1=1,3 
145 XCA(l+l) = X2GA(I) 

XCA(5) = ACM 
DO 150 1=1,3 
150 XCA(I45) = Y1GA(I) 

XCA(9 ) = ARP 
XDA(l) = ALP 
DO 155 1=1,3 
155 XDA(I+1) = X1GA(I) 

XDA(5) = ACM 
DO 160 1=1,3 
160 XDA( 1+5 ) = Y1GA(I) 

XDA(9) = ARP 
GO TO 500 

200 IF (S(2) .LT.0.0) GO TO 300 
XAA(l) = ALP 
DO 225 1=1,3 
225 XAA(I+1) = X2GA(I) 

XAA( 5 ) = ACM 
DO 230 1=1,3 
230 XAA( 1+5 ) = YlGA(l) 

XAA(9) = ARP 

XBA(l) = ALP 

DO 235 1=1,3 
235 XBA(I+1) = X2GA( I ) 

XBA( 5 ) = ACM 
DO 240 1=1,3 
240 XBA( 1+5 ) = Y2GA( I ) 

XBA{9 ) = ARP 

XCA(l) = ALP 

DO 245 1=1,3 

245 XCA( 1+1 ) = X1GA(I) 

XCA( 5) = ACM 
DO 250 1=1,3 
250 XCA( 1+5 ) = Y1GA{I) 

XCA(9 ) = ARP 

XDA(l) = ALP 

DO 255 1=1,3 
255 XDA( 1+1 ) = X1GA(I) 

XDA(5) = ACM 
DO 260 1=1,3 
260 XDA( 1+5 ) = Y2GA(I) 

XDA(9 ) = ARP 
GO TO 500 
300 XAA(l) = ALP 
DO 325 1=1,3 
325 XAA( 1+1 ) = XlGA{ I ) 

XAA(5) = ACM 
DO 330 1=1,3 
330 XAA( 1+5 ) = Y2GA{I) 

XAA( 9 ) = ARP 
XBA(l) = ALP 
DO 335 1=1,3 
335 XRA( 1+1 ) = X1GA(I) 

XBA{5) = ACM 
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IT) 340 1=1,3 

340 XBA( 1+5 ) = YlGA(l) 

XBA(9) = ARP 
XCA(l) = ALP 
DO 345 1=1,3 
345 XGMl+1) = X7GA( I ) 

XCA(5) = ACM 
DO 350 1=1,3 
350 XCMl+5) = Y?GA( I ) 

XCA (9 ) = ARP 
XDA(l) - ALP 
DO 355 1=1,3 
355 XDA( 1+1 ) = X2GA( I ) 

XDA(5 ) = ACM 
DO 300 1=1,3 

360 XDA( 1+5 ) = YlGA(l) 

XDA{9 ) = ARP 
C TYPE 905, XDA 

500 CALL SYMBOL (XB-. 5, YB-. 25, .075, XBA , 0. 0, 17 ) 
CALL SYMBOL (XA-. 75, YA+. 25, . 075,XAA, 0. 0, 17) 
CALL SYMBOL (XC-. 75, YC+. 25, . 075, XCA, 0. 0, 1 7 ) 
CALL SYMBOL (XD-.75, YD-.25, .075, XDA, 0.0, 17) 
CALL NLMCER (XO, YKAX+. 1 , . 1 , ZCKAX, 0 . 0, J ) 

CALL PL0T(0. , 0. , 499° ) 

.STOP 

END 
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APPENDIX B 

TYPICAL PLOTS OF FLUX-MAPPER DATA 
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Figure B-l. A Typical Flux 
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Figure B-2. A Typical Contour Map at a Concentrator Focal Plane 
(Contour software not listed in this report, but 
available from SNLA.) 
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